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CHAPTER  I 


INTRODUCTION 


1.1  Objective 

An  approximate  expression  for  the  far  zone  field  scattered  by  the  vertex  of 
a  finite  perfectly  conducting  wedge  is  obtained.  The  solution  is  cast  in  the  form 
of  the  UTD  and  is  based  on  asymptotic  eqmvalent  currents  found  using  modified 
PTD  concepts.  The  faces  of  the  wedge  must  be  flat  (the  normal  to  each  individual 
face  is  a  constant  everywhere  on  the  face  except  at  the  edge)  and  the  edges  must 
be  straight.  For  plane  wave  incidence  from  an  arbitrary  direction  the  first  order 
contribution  from  each  vertex  to  the  far  zone  scattered  field  is  obtained. 

1.2  Motivation 

Since  difiraction  is  a  local  phenomena  at  high  frequencies  the  results  obtained 
for  a  finite  wedge  may  be  applied  to  much  more  complex  bodies  made  up  of  simple 
shapes.  The  field  scattered  by  a  three-dimensional  shape  constructed  from  flat 
plates  may  be  approximated  to  first  order  u  the  sum  of  the  contributions  from 
each  individual  corner.  The  first  order  solution  should  be  reasonably  accurate  in 
or  near  the  specular  regions  as  long  as  the  object  is  convex.  A  convex  body  is 
defined  here  as  a  dosed  surface  made  up  of  fiat  plates  such  that  all  of  the  exterior 
wedge  angles,  taken  between  faces  and  exterior  to  the  surface,  are  greater  than  180 
degrees.  A  simple  example  of  an  object  that  does  not  meet  this  requirement  is  a 
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corner  reflector.  The  effect  of  shadowing  of  the  faces  by  other  parts  of  the  object, 
however,  must  be  taken  into  account.  A  brief  description  of  how  this  is  done  is 
given  in  Section  5.2. 

1.3  Problem  History 

There  are  many  approximate  solutions  to  the  scattered  field  from  a  finite 
perfectly  conducting  wedge.  Some  of  these  include  the  Physical  Optics  approx¬ 
imation,  the  Method  of  Equivalent  Currents  [1],  the  previous  Comer  Diffraction 
coefficients  [2],  the  Uniform  Geometrical  Theory  of  Diffraction  [3], and  the  Physical 
Theory  of  Di&action  [4].  The  Physical  Optics,  Method  of  Equivalent  Currents, 
rnd  Physical  'THeoi^  of  Diffraction  will  be  described  briefly  in  the  next  chap¬ 
ter.  The  previous  Comer  Diffraction  Coefficient  was  heuristically  derived  from  the 
equivalent  run  mts  of  Ryan  and  Peters  [1]  and  the  Uniform  Geometrical  Theory 
of  Diffraction  [3].  Although  it  ^ves  good  results  in  the  backscatter  region,  it  has 
been  found  to  (pve  unsatisfactory)’  results  for  certain  bistatic  cases.  Specifically,  it 
has  problems  at  the  so-called  false  shadow  boundaries  which  will  be  illustrated  in 
Chapter  II. 

1.4  Outline  of  Solution 

It  is  assumed  that  the  incident  field,  and  therefore  the  t::attered  field,  is  a 
time  harmonic  field  with  time  dependence  ’riven  by  which  will  be  suppressed 
throughout  this  report.  The  solution  will  be  ^  ased  on  the  PTD  and  cast  into  the 
form  of  the  method  of  equivalent  currmts  and  then  into  diffraction  coefficients. 
This  will  he  done  a£  follows.  The  actual  currents  on  each  face  of  the  plate  will  be 
approximated  by  the  PTD  currents.  The  resulting  double  integral  ovei  the  surface 
Is  then  reduced  to  a  line  integral  along  the  edge  by  doing  the  integration  over 
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one  coordinate  in  closed  form.  This  is  done  by  taking  the  asymptotic  endpoint 
contribution  in  one  case  and  by  the  application  of  some  algebraic  manipulation 
followed  by  the  use  of  Stoke’s  theorem  in  the  other  case.  The  remaining  integral 
«  along  the  edge  is  written  in  the  form  of  a  radiation  integral  so  that  the  equivalent 

currents  may  be  identified.  Then  this  integral  is  evaluated  using  the  method  of 

f 

stationsLTy  phase  to  obtain  the  contribution  &om  each  comer. 

1.6  Notation,  Abbreviations  and  Symbols 

This  section  contains  some  information  on  the  notation,  abbreviations,  and 
symbols  that  are  used  in  this  report.  Normally  the  term  equivalent  currents  refers 
to  the  surface  currents  used  in  the  equivalence  theorem.  In  this  report,  however,  the 
terms  equivalent  edge  currents  and  equivalent  currents  will  be  used  interchangeably 
to  refer  to  equivalent  edge  currents. 

1.6.1  Symbols 

These  are  some  of  the  symbols  used.  Most  of  them  are  also  defined  as  they 
first  appear  in  the  text. 

0  =  diffracted  ray  elevation  angle  in  edge  fixed  coordinates 

A 

13  —  diffracted  ray  elevation  unit  vector  in  edge  fixed  coordinates 
0'  =  incident  ray  elevation  angle  in  edge  fixed  coordinates 

-  incident  ray  elevation  unit  vector  in  edge  fixed  coordinates 
A  =  free  space  wavelength  of  the  time  harmonic  fields 
•  ^  =  diffracted  ray  azimuth  angle  in  edge  fixed  coordinates 

or  in  pattern  coordinates  (depending  on  context) 

^  —  diffracted  ray  azimuth  unit  vector  in  edge  fixed 
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coordinates  or  in  pattern  coordinates  (depending  on  context) 
=  incident  ray  azimuth  angle  in  edge  fixed  coordinates 
or  in  pattern  coordinates  (depending  on  context) 

=  incident  ray  azimuth  unit  vector  in  edge  fixed 

coordinates  or  in  pattern  coordinates  (depending  on  context) 
=  Radar  Cross  Section 

=  diffracted  ray  elevation  angle  in  pattern  coordinates 
=  diffracted  ray  elevation  unit  vector  in  pattern  coordinates 
=  incident  ray  elevation  angle  in  pattern  coordinates 
=  incident  ray  elevation  vector  in  pattern  coordinates 
=  angular  frequency 
=  bi-normal  unit  vector 
=  soft,  hard  diffraction  coefficients 
=  soft,  hard  corner  diffraction  coefficients 
=  cross  polarized  corner  diffraction  coefficient 
(hard  incidence,  soft  diffracted) 

=  corner  diffracted  field 
=  diffracted  field 
=  incident  field 

=  components  of  the  incident  field  tangent  to  the  edge 
=  transition  function 


i  = 


le,m  =  electric  equivalent  currents 


$ 
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k  =  27r/A 

M  =  magnetic  eqmvalent  current 
n  =  wedge  angle  parameter 
n  =  unit  vector  normal  to  the  face  of  the  wedge 
<  =  unit  vector  tangent  to  the  edge 
s'  =  incident  ray  direction 
s  =  diffracted  ray  direction 
Kq  =  admitt2mce  of  free  space 
Zq  =  impedance  of  free  space 

Components  of  the  equivalent  currents  will  be  designated  by  placing  appropriate 
superscripts  on  the  symbols  and  M. 

1.5.2  Abbreviations 


GO  =  Geometrical  Optics 
GTD  =  Geometrical  Theory  of  Diffraction 
MM  =  Moment  Method 
PEC  =  Perfect  Electrical  Conductor 
PO  =  Physical  Optics 
PTD  =  Physical  Theory  of  Diffraction 
RCS  =  Radar  Cross  Section 
UTD  =  Uniform  Geometricid  Theory  of  Diffraction 
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1.5.3  Edge  Fixed  Coordinate  System 

An  edge  fixed  coordinate  system  will  be  used  throughout  the  discussion  of  the 
equivalent  currents  and  the  comer  diffraction  coefficient.  As  shown  in  Figure  1, 
the  origin  of  the  right  hand  system  (b,n,t)  is  placed  on  the  edge.  The  point  on 
the  edge  where  the  oripn  is  placed  will  depend  on  the  problem  being  considered 
and  will  be  specified.  The  n  direction  is  chosen  so  n  is  normal  to  the  face  under 
consideration  and  directed  outward  &om  the  interior  of  the  body.  In  the  case  of  a 
flot  plate  the  outward  normal  is  ambiguous,  and  the  positive  n  direction  may  be 
arbitrarily  chosen  in  either  of  the  two  directions.  The  (  axis  is  chosen  tangent  to 
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the  edge.  The  positive  i  direction  is  chosen  such  that  the  positive  b  axis  {b  =  hxi) 
will  lie  on  the  wedge  face  under  consideration.  Since  only  one  face  is  considered  at 
a  time  this  results  in  a  convenient  and  unambiguous  coordinate  system  for  all  of 
the  problems  considered  here. 

The  incident  ray  is  described  by  the  spherical  coordinate  system  {s' ,13' ,<!>') 
based  on  an  axis  in  the  ~t  direction.  The  incident  field  may,  therefore,  be  described 

A  A 

as  a  sum  of  its  components  in  the  and  <l>'  directions  since  it  is  always  assumed  to 
be  a  plane  wave  in  this  report.  Similarly  the  diffracted  rays  are  described  by  the 
spherical  coordinate  system  {s,P,<j))  based  on  an  axis  in  the  positive  i  direction. 
The  radial  component  of  the  diffracted  field  will  be  zero  in  the  far  zone,  resulting 
in  a  field  which  may  be  expressed  in  terms  of  its  and  ^  components.  The  angles 
and  unit  vectors  are  shown  in  Figure  1. 

1.5.4  Radar  Cross  Section 


The  numerical  results  are  expressed  in  terms  of  Radar  Cross  Section  through¬ 
out  this  report.  Radar  Gross  Section  is  the  mathematical  area  that  intercepts  the 
power  from  the  incident  wave  on  the  target,  which  if  scattered  in  an  isotropic  man¬ 
ner  would  produce  the  same  scattered  power  density  as  produced  by  the  actual 
target  assuming  that  both  the  transmitter  and  the  receiver  are  in  the  far  zone  of 
the  target;  i.e.  ri,  ^2  >  2D^JX,  where  D  is  the  maximum  dimension  of  the  target 
and  ri,r2  arc  defined  in  Figure  2.  Letting  5^=  power  density  incident  on  the 
target  from  the  transmitter  (constant  over  the  target  for  plane  wave  incidence); 

power  density  scattered  from  the  target  toward  the  receiver  (constant,  see 
S*)  results  in 


Un. 

’•2-00  iirrl 


(1.1) 
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Figure  2:  Definition  of  Distances  rj  and  r2. 


Solving  for  <r  gives 


Hm 


(1.2) 


Since  for  most  problems  of  practical  importance  and  all  of  the  problems  done  here 
the  medium  is  free  space  the  power  densities  are  simply  related  to  the  fields  by 

IBI* 


S  =  ^  =  \H\^Z^ 


(1.3) 


where  |£l|  and  |if|  are  the  rms  values  of  the  electric  and  magnetic  fields,  respec¬ 
tively,  at  the  point  of  interest,  and  Zq  is  the  impedance  of  free  space.  Therefore 
the  RCS  reduces  to 


<r  =  lim  4irri — =  lim  4irri 

r2-*oo  ^  I  _--l2  r«-»oo  * 


l^fci2 


£?* 


r2-»oo 


if* 


ir 


(1.4) 


where  and  |ir*®|  are  the  magnitude  of  the  scattered  electric  and  magnetic 


fields,  respectively,  at  the  receiver  and  and  |jl7*|  are  the  magnitude  of  the 
incident  electric  and  magnetic  fields  at  the  target. 
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1.6  Ptittern  Coordinate  System 


The  usual  spherical  coordinate  system  is  used  for  all  of  the  results.  The  source 
is  located  in  the  far  zone  direction  specified  by  the  spherical  angles  and  (or 
sometimes  by  and  ^*).  Similarly  the  observation  point  is  located  in  the  far  zone 
direction  specified  by  the  sphericsd  angles  6  and  <ff.  In  the  special  case  of  backscatter 
both  the  source  and  receiver  locations  will  be  indicated  by  6  and  Notice  that  ^ 
and  (p'  are  also  used  in  defining  the  observation  and  source  directions,  respectively, 
in  edge  fixed  coordinates,  but  it  should  be  obvious  from  the  context  which  angle  is 
being  identified.  The  RCS  is  given  in  terms  of  the  polarization  of  the  transmitted 
and  received  fields.  The  following  definitions  will  be  used  throughout  this  report 


<Ts9  =  4irr2 


or^  =  4irr| 
=  4wf  f 


-  4jrr2 


2 

4 

f 

4 

1 

2 

4 

5 

2 

I 

(1.5) 


(1.6) 

(1.7) 


(1.8) 


where 


and  |£^j  are  the  magnitude  of  the  4>'  and  9'  components,  respectively, 
of  the  incident  field  at  the  target.  Sinularly  and  are  the  ^  and  B 
components  of  the  scattered  field  at  the  receiver.  The  first  subscript  on  (r  refers 
to  the  polarization  of  the  received  field;  whereas,  the  second  one  refers  to  the 
polarization  of  the  incident  field,  both  in  pattern  coordinates.  For  the  general  case 
of  bistatic  scattering,  the  notation  will  remiun  the  same  even  though  9  and  9'  (and 
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^  Old  “  different  directions.  All  results  shown  nre  given  in  terms  of  either 

dB  relative  to  a  aqnare  meter  or  dB  relative  to  a  square  wavelength  depending  on 
the  specific  problem.  They  are  simply  related  to  the  above  expresrion.  by 

<r(mdB)  =  101ogio<r 


or 


<r(m  dB)  ~  20logio 


+  10logio(47r) 


for  each  of  the  combinations  of  transmitting  and  receiving  polarizations. 


(1.10) 
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CHAPTER  n 


THEORETICAL  BACKGROUND 


2.1  Introduction 

The  high  frequency  methods  used  to  derive  the  new  far  sone  corner  diffrac¬ 
tion  coefficients  are  described  briefly  here.  Geometrical  Optics  (GO)  and  Physical 
Optics  (PO)  are  discussed  first  since  both  of  these  are  important  in  the  implemen¬ 
tation  of  the  new  solution.  The  basic  concepts  of  the  Physical  Theory  of  Diffraction 
(PTD)  are  then  given.  Following  this  is  a  description  of  the  differences  between 
the  Method  of  Equivalent  Currents  as  proposed  by  Ryan  and  Peters  [1]  and  the 
method  used  to  derive  the  new  equivalent  currents.  A  corner  diffraction  coefficient 
based  on  the  old  equivalent  currents  is  also  given.  The  results  from  the  method  of 
stationary  phase  are  then  described  briefly. 

2.2  Geometrical  Optics 

A  brief  summary  of  Geometrical  Optics  (GO)  is  given  here.  Further  informa¬ 
tion  on  the  subject  may  be  found  in  references  (5]  -  [8],  and  the  basic  principles 
are  covered  in  many  antennas  textbooks,  such  as  [9].  The  propagation  of  elec¬ 
tromagnetic  energy  through  isotrofuc,  lossless  media  can  be  described  using  GO. 
It  has  long  been  known  that  at  high  frequencies  electromagnetic  energy  can  be 
viewed  as  traveling  along  well  defined  paths  known  as  rays.  The  ray  paths,  in 
any  continuous  medium,  may  be  determined  using  Fermat’s  principle  which  states 
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the  energy  will  flow  along  the  path  of  shortest  electrical  length  between  any  two 
points  in  the  medium.  The  shortest  electrical  path  is  the  path  which  results  in  the 
shortest  propagation  time  between  the  two  points.  (Hirthermore,  these  rays  are 
orthogonal  to  surfaces  of  constant  phase.  If  the  media  is  also  homogeneous,  which 
is  true  in  this  case,  the  ray  paths  are  straight  lines.  Assuming  the  fleld  is  known  at 
some  constant  phase  reference  surface,  it  can  be  calculated  at  any  point  away  from 
caustics  using  conservation  of  energy.  Referring  to  Figure  3,  if  the  fleld  is  known 
at  some  constant  phase  surface  HAq  along  with  the  principle  radii  of  curvature,  p\ 
and  P2»  of  the  surface  dAo  then  the  field  at  dA  may  be  calculated  using  the  conser> 
vation  of  energy  and  the  above  assumption  that  the  electromagnetic  energy  travels 
in  straight  lines.  Since  the  energy  flux  at  both  of  the  surfaces  is  proportional  to 
the  square  iA  the  field,  it  follows  from  the  conservation  of  energy  that 

|jg(0)|*«o  =  ISwf'IA . 

FVoro  gemnetric  comnderations  it  can  be  shown  that 

^  0(p2  0 

dAo  FlPl 


(2.1) 

(2.2) 
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Combining  these  two  equations  gives 


i(<)l  =  |e(o)| 


P1P2 

(Pl  +  ^)(P2  + 

(2.3) 


where  it  is  assumed  that  the  positive  t  direction  is  in  the  direction  of  propagation 
and  the  principle  radii  of  curvature  p\  and  p2  may  be  both  positive  (concave  wave 
front),  both  negative  (convex  wave  front),  or  opposite  in  sign  (saddle  wave  front). 
The  complete  expression  for  the  field,  one  that  includes  the  phase  as  well  as  the 
magnitude  of  the  field,  must  be  obtained  from  the  asymptotic  solution  of  Maxwell’s 
equations  and  is  given  v 


ii'(^)  -=  E{Q) 


1 


PiP2 


(Pl  +  ^)(P2  +  ^) 


-i(W-mf) 


(2.4) 


where  m  is  the  number  of  caustics  the  ray  passes  through  going  from  0  to  £.  For 
example,  ra=0  for  points  to  the  right  of  the  caustic  labeled  3-4  as  shown  in  Figure  3, 
m=l  for  pcims  between  1-2  and  3-4,  and  m=:?  for  points  to  the  left  of  1-2.  It  has 
been  assumed  that  the  wave  is  traveling  from  left  to  right  so  i  is  positive  to  the 
right  of  0  and  negative  to  the  left  of  O.  GO  obviously  fails  at  the  caustics  where 
it  predicts  that  the  field  becomes  infinite. 

Only  two  special  cases  are  of  interest  here.  Letting  pi,P2  -*  oo  and  replacing 
£  with  8  gjves;  £(«)  =  £(0)e~^^*  which  is  the  well  known  result  for  a  plane  wave 
propagating  in  the  positive  s  direction  where  E(0)  is  the  field  at  any  convenient 
constant  phase  plane,  and  s  is  the  distance  from  the  reference  plane  to  the  point 
of  interest.  The  spherical  wave  is  the  second  case  of  interest  here.  In  this  case, 


P1=P2-P  giving 


^(s)  =  E(0) 


p  +  s 


letting  the  reference  point  approach  the  caustic  and  assuming  pE(0) 


(2.5) 
As  as 
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n 


Figure  4:  Projection  of  rays  on  to  b-n  plane. 


p  gives 


E{a)  =  As- - 


(2.6) 


which  is  the  familiar  spherical  wave  result.  Notice  that  m=0  in  both  cases  since 
no  caustics  are  crossed. 

2.S  Geometrical  Theory  of  Diffraction 


If  Fermat’s  principle  is  extended  to  include  refraction,  reflection,  and  diffrac¬ 
tion  points  then  GO  may  be  extended  to  media  with  discontinuities  in  electrical 
characteristics.  This  means  that  not  only  is  the  held  from  the  direct  or  line-of- 
sight  path  included,  but  also  paths  which  include  points  on  the  surface  where  the 
media  is  discontinuous.  In  this  case  the  discontinuity  is  a  wedge  with  perfectly 
conducting  flat  faces.  The  edge-fixed  coordinate  system  described  earlier  is  used. 
The  projection  of  the  rays  onto  the  b-n  plane  is  shown  in  Figure  4.  For  a  fixed 
source,  the  observation  point  may  be  in  one  of  three  regions  around  the  wedge.  In 
region  1  (0  <  ^  <  ir  -  ^'),  all  Uiree  rays  contribute  to  the  total  field.  In  region 
11  (x  -  <  ^  <  X  +  only  the  direct  and  diffracted  fields  contribute  to  the 
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OBSERVATION 

POINT 


Figure  S:  Keller’s  cone  of  diffracted  rays. 

total  field  since  no  point  on  the  wedge  satisfies  the  law  of  reflection.  In  region  III 
(w  4-  <  <f>  <  tmt),  only  the  diffracted  field  contributes  to  the  total  field.  No¬ 

tice  that,  in  general,  the  reflected  field  is  discontinuous  at  the  reflection  shadow 
boundary  (^-f  =  ir),  and  the  direct  fi*H  is  discontinuous  at  the  incident  shadow 

boundary  -  <f>'  =  w).  Keller  [10]  derived  an  expression  for  the  diffracted  field 
in  terms  of  the  field  incident  on  the  edge  and  showed  that  the  major  contribution 
comes  hrom  a  ;dngle  point  (the  diffraction  point).  The  condition  0  ~  0*  ib  met  at 
the  diffraction  point.  The  cone  of  rays  that  satisfy  this  condition  for  a  given  point 
along  .  ue  edge  is  known  as  the  KeUe>  cone  (see  Figure  5).  Keller’s  results  for  the 
diffracted  field  were  only  for  observation  points  away  from  the  shadow  boundaries. 
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Kouyounyian  and  Pathak  (3j  later  formulated  the  Uniform  Geometrical  Theory  of 
Dilfraction  (UTD)  which  extended  the  solution  to  all  regions  of  space.  The  results, 
in  terms  of  the  edge  fixed  coordinates  shown  earlier  in  Figure  1,  for  a  plane  vrave 
incident  from  the  direction  and  a  far  field  observation  direction  at  (/3,(^)  are 

given  by 

r  T  I  II  I 


Ds 

0 

y*. 

0 

Dh 

4(<3.) 

(2.7) 


where 


^s,h  ~ 


-e-n 


+  («>  -  »')M 

2n 


Qt  is  the  point  on  the  wedge  which  satisfies  the  condition  0  =  0',  Ep,(Qe}  is 
the  component  of  the  incident  field  in  the  0'  direction  at  Qe  i® 

component  of  the  incident  field  in  the  direction  at  Qg  (see  Figure  5).  The 
diffracted  field  is  described  in  terms  of  transverse  components,  0  and  6,  at  a  far 
sone  distance  s.  The  first  and  third  cotangent  terms  have  been  associated  with 
the  0-face  incident  and  reflection  boundaries,  respectively.  Similarly  the  second 
and  fourth  cotangent  terms  are  associated  with  the  N-face  incident  and  reflection 
boundaries. 

2.4  Physical  Optics 

Physical  Optics  is  a  widely  used  technique  for  finding  the  field  scattered  by  a 
perfectly  conducting  body  of  arbitrary  shape.  The  physical  optics  approximation 
to  the  surface  current  is  given  in  [11]  by; 

2nx  S'  \  in  Hi  regions 


JpO  = 


(2.8) 


0  ;  in  shadow  regions 
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Figure  6:  PO  Currents  on  an  arbitrary  scatterer. 
where  the  illuminated  and  shadowed  regions  are  determined  using  ray  optics,  and 

maf  » 

is  the  incident  magnetic  field  on  the  surface  of  the  perfectly  conducting  body 
and  is  approximated  by  the  GO  incident  held  in  this  case.  Since  only  far  zone 
sources  are  considered  here  the  incident  held  will  always  be  locally  plane  and  the 
rays  will  be  parallel.  A  simple  2'D  example  illustrating  the  lit  and  shadowed 
regions  is  given  in  Figure  6. 

The  scattered  held  is  found  using  the  far  zone  radiation  integral  which  is 
dehned  by 

S'(>)  =  /j,  .  (2.9) 

Plugging  in  the  PO  approximation  to  the  surface  current  gives: 

g»(,)  =  2n  X  (2.10) 
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where  S'  is  the  illuminated  region  on  the  scattering  surface  and  and  r'  are 
defined  in  Figure  6. 

2.5  Physical  Theory  of  Diffraction 

The  Physical  Theory  of  Diffraction  was  first  widely  used  for  scattering  prob¬ 
lems  by  Ufimtsev  [4].  Additional  currents  are  added  to  the  Physical  Optics  currents 
so  that  the  solution  is  valid  for  larger  regions  of  space  and  for  smaller  scatter- 
ers  than  is  possible  with  either  GO  or  PO.  Letting  the  total  surface  current  be 
J  —  Jpo  Jf  where  is  the  correction  to  the  PO  current  which  when  added  to 
the  PO  current  would  approximate  the  actual  current  on  the  scatterer.  Since  the 
PO  current  is  derived  from  the  assumption  that  the  surface  may  be  locally  approx¬ 
imated  as  an  infinite  plane  tangent  to  the  surface,  it  is  expected  that  Jf  could  be 
described  as  the  current  due  to  the  deviation  of  the  smfacc  from  an  infinite  plane. 
The  sharp  edge  or  wedge  is  the  only  deviation  from  a  planar  surface  considered 
here.  Ufimtsev  speculated  that  this  correction  to  the  PO  current  would  only  be 
important  near  the  edges  of  the  scatterer  in  this  case. 

Assuming  plane  wave  incidence,  the  magnitude  of  the  PO  component  of  the 
current  is  constant  on  the  illuminated  region  of  a  planar  scatterer  so  Ufimtsev 
named  this  the  uniform  component  of  the  current.  The  magnitude  of  the  cor¬ 
rection  current  is  obviously  not  constant  on  the  surface  and  is  thus  called  the 
nonuniform  component  of  the  current.  In  the  case  of  a  wedge,  Ufimtsev  refers  to 
the  field  radiated  by  the  nonuniform  part  of  the  current  as  an  elementary  edge 
wave  or  rimply  an  edge  wave  since  it  is  produced  by  the  nonuniform  part  of  the 
current  which  is  concentrated  along  the  edge  and  rapidly  decreases  away  from  the 
edge.  In  his  latest  paper  [12],  however,  it  seems  as  though  he  refers  to  the  scattered 
field  as  the  total  edge  wave  and  the  field  due  to  the  correction  current,  Jy,  on  an 
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infinitesimal  strip  of  the  wedge  as  the  elementary  edge  wave.  In  this  case  the  only 
contribution  to  the  nonuniform  current  Jf  h  due  to  the  edge  discontinuities.  If 
the  edges  are  long  in  terms  of  the  wavelength  and  have  radii  of  curvature  which 
are  also  long  in  terms  of  the  wavelength  the  total  current  may  be  approximated 
by  the  current  on  an  infinite  wedge  tangent  to  the  scatterer  at  the  point  of  inter¬ 
est.  The  nonuniform  component  of  the  current  may  then  be  found  by  subtracting 
the  GO  current  from  this  approximation  to  the  total  current.  Ufimtsev  calls  the 
approximation  to  the  nonuniform  current  given  by  the  infinite  wedge  the  fringe 
current,  and  the  field  produced  by  this  current  the  fringing  field.  Ufimtsev  never 
actually  found  the  fringe  current,  instead  he  found  the  fringing  field  by  indirect 
considerations.  In  recent  years  both  Michaeli  [13,14]  and  Ufimtsev  [12,15]  have  de¬ 
rived  expressions  for  the  fringing  field  involving  expressions  for  the  fringe  current 
on  the  wedge.  Michaeli’s  expressions  for  the  fringe  current  will  be  used  in  the  later 
development  of  the  corner  diffraction  coefficient. 

2.6  Method  of  Equivalent  Currents 

The  method  of  equivalent  currents  was  originally  used  by  Millar  in  [16]  -  [18]. 
Ryan  and  Peters  [1]  used  the  method  of  equivalent  currents  to  find  the  diffracted 
field  in  the  region  of  caustics  of  the  UTD  diffracted  field.  Ryan  and  Peters  equiva¬ 
lent  currents  have  been  used  to  find  the  original  comer  diffraction  solution  shown 
in  the  next  section.  Ryan  and  Peters  equivalent  currents  are  discussed  here  since 
they  may  be  used  to  illustrate  some  important  points  about  the  use  of  equiva¬ 
lent  currents  in  general.  A  brief  description  of  the  methods  of  finding  the  new 
equivalent  currents  used  in  the  corner  diffraction  coefficient  will  be  given. 

Ryan  and  Peters  compared  the  asymptotic  approximation  to  the  far  zone 
diffracted  fields  from  an  infmite  wedge  illuminated  by  a  plane  wave  to  the  expres- 
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sions  for  the  far  zone  fields  radiated  by  infinite  magnetic  and  electric  line  sources. 
Equating  the  two  fields,  they  obtsdned  expressions  for  the  equivalent  edge  currents 
which  depend  on  the  direction  of  incidence  and  the  direction  of  observation.  The 
field  scattered  by  an  arbitrary  three-dimensional  body  made  up  of  smoothly  curved 
faces  terminated  in  sharp  edges  is  found  by  numerically  evaluating  the  usual  radi¬ 
ation  integral  along  the  edges  of  the  scatterer  using  the  equivalent  currents.  These 
equivalent  edge  currents  are  obviously  not  physical  currents  since  they  are  func¬ 
tions  of  the  observation  point.  They  are  a  means  of  mathematically  simplifying 
the  problem  of  finding  the  diffracted  field. 

It  was  shown  later  by  Northrop  [19]  and  Sikta,  et  ai  [2]  that  if  the  surfaces 
were  represented  by  very  thin  strips  oriented  in  the  proper  direction,  the  solution 
was  more  accurate  for  three-dimensional  problems.  The  strips  are  oriented  in  a 
direction  such  that  the  angle  between  the  vector  tangent  to  the  edge  (t)  and  the 
strip  is  given  by 

iaa9$  = cos<l>'  0  <  <  tt  (2.11) 


as  shown  in  Figure  7.  This  orients  the  strips  so  that  the  incident  ray  is  normal  to 
the  infinitesimal  edge  on  each  strip. 

Later,  it  was  found  by  Marhefka  [20]  that  for  the  case  of  bistatic  scatter¬ 
ing  changing  the  orientation  of  the  strips,  in  general  improved  the  results.  The 
modification  used  is  simply  to  replace  e*  given  in  [2]  by: 


-4 

e*  = 


X 

1 

f  (I-i)xn 

A 

1 

X 

[  (/  -  s)  X  n 

®P 

1  < 

(2.12) 


where  s  is  the  diffraction  direction,  /  is  the  incident  ray  direction,  n  is  the  normal 
to  the  face,  and  ip  is  the  edge  vector.  This  expression  reduces  to  the  expression  for 
it  given  in  [2]  for  the  special  case  of  backscatter.  This  replacement  is  equivalent 


20 


Figure  7;  Strips  used  in  Ryan  and  Peters  Equivalent  Currents  by  Sikta  [2j. 
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to  placing  the  strips  such  that: 


cot  fig 


cosjS  —  cos  10^ 


sin  cos  4>  +  sin  cos 


0  <  <  TT 


(2.13) 


where  Bg  is  measured  as  illustrated  in  Figure  7.  It  is  interesting  to  note  that  this 
is  the  strip  orientation  associated  with  the  PO  surface  integral  (see  Section  4.7). 

As  a  result  of  the  two-dimensional  nature  of  the  solution  and  the  orientation  of 
the  strips,  the  scattered  field  predicted  by  the  equivalent  currents  is  singular  at  the 
so-called  false  shadow  boundaries.  The  false  shadow  boundaries  are  the  regions 
of  space  where  the  two-dimensional  problem  goes  through  a  shadow  boundary 
~  Tr)»  but  the  three-dimensional  problem  does  not  ^  (3').  The  major 
reason  for  this  problem  is  the  two-dimensional  nature  of  the  solution.  A  simple 
example  is  given  next  illustrating  this  property. 

The  bistatic  scattering  &om  a  flat  plate  two  wavelengths  on  a  side  is  calculated 
using  the  equivalent  currents  given  in  (2)  with  the  replacement  of  e*  as  described 
above.  A  fixed  source  is  placed  at  ^  =  45®  and  =  0®  as  shown  in  Figure  8. 
The  pattern  is  taken  in  the  ^  =  60®  plane  and  the  results  are  illustrated  for  the 
entire  scattering  matrix.  The  results  using  the  previous  equivalent  currents  are 
compared  with  Method  of  Moments  results  in  Figures  9  to  12  for  the  four  different 
combinations  of  incident  field  and  diffracted  field  polarizations.  In  all  four  cases, 
the  previous  equivalent  current  solution  is  singular  at  the  false  shadow  boundaries 
{6  Si  240®  and  6  Si  300®  in  this  case).  The  spikes  are  due  to  singularities  in 
the  contributions  from  both  the  front  and  back  edges  which  have  false  shadow 
boundaries  at  both  9  »  240®  and  9  ss  300®. 

The  idea  of  using  infinitesimal  strips  to  represent  a  fiat  face  is  used  again  in 
deriving  the  new  equivalent  currents;  however,  the  strips  are  oriented  in  a  different 
direction  in  finding  the  new  equivalent  currents.  A  different  approach  is  used  to 
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Figure  8:  FUse  Shadow  Boundary  example  geometry  (2A  square  plate  in  the  x-y 

plane). 
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(dB/  X’) 


0  (deg.) 


Figure  9:  Co-polarired  RCS  in  the  ^  =  60°  plane  of  a  2A  square  plate  with  a  6* 
polarized  fixed  source  at  C*  =  45°,  =0°. 


Figure  10:  Co-polarized  RCS  in  the  =  00°  plane  of  a  2A  square  plate  with  a  ^ 
polarized  fixed  source  at  9*  =  45°,  =  0°. 
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0  (deg.) 

Figure  11;  Gross-polarized  RCS  in  the  ^  =  60®  plane  of  a  2 A  square  plate  with  a 
¥  polarized  fixed  source  at  =  45®,  =  0®. 


Figure  12:  Cross- polarized  RCS  in  the  4>  ~  plane  of  a  2A  square  plate  with  a 
polarized  fixed  source  at  =  45®,  =0®. 
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find  the  new  equivalent  edge  currents,  and  different  expressions  are  obtained.  The 
scattered  field  is  found  from  the  new  equivalent  edge  currents  in  the  same  way  as  it 
is  found  from  Ryan  and  Peters  previous  expressions.  The  new  equivalent  currents, 
however,  predict  smooth  fields  at  the  false  shadow  boundaries.  This  is  illustrated  in 
Section  5.3.7  where  the  above  example  is  repeated  for  the  new  equivalent  currents. 

2.7  Previous  Corner  Diffraction  Coefficient 


A  diffraction  coefficient  for  a  cora'^r  formed  by  the  intersection  of  two  straight 
edges  was  derived  by  Burnside  and  Pathak  [2].  It  is  based  on  the  asymptotic 
evaluation  of  the  radiation  integral  containing  the  equivalent  currents  of  Ryan 
and  Peters  [1]  The  result  was  then  empirically  modified  so  that  the  diffraction 
coefficient  would  not  change  sign  abruptly  as  it  passes  through  the  false  shadow 
boundaries.  It  was  derived  for  spherical  wave  incidence  and  remains  valid  for  cases 
when  the  diffraction  point  is  near  the  corner  since  the  integral  was  evaluated  for 
a  saddle  point  near  an  end  point;  however,  only  the  far  zone  result  is  shown  here. 
The  corner  diffracted  field  due  to  one  corner  and  one  edge  in  the  case  of  plane 
wave  incidence  and  a  far  zone  receiver  is  given  by 
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a(/9)  =  2co6^  =  2008^ 


where  is  the  integer  which  most  nearly  satisfies  2nirN^  —  ^  =  :f  tt,  and 


(2.16) 


where  the  angles  are  shown  in  Figure  13.  The  sign  on  the  diffraction  coefficient 
may  be  plus  or  minus  depending  on  which  endpoint  of  the  edge  is  being  considered. 
Choosing  the  correct  sign  is  discussed  in  Section  3.5. 

The  heuristic  factor  added  to  the  corner  diffraction  coefficients  changes  the 
behavior  of  the  solution  at  the  false  shadow  boundaries,  as  illustrated  in  the  next 
'  example.  The  bistatic  scattering  e;:ample  from  Section  2.6  (2A  square  plate,  source 
fixed  at  0*  —  45°  and  =  0°,  and  ^  =  60°  pattern)  is  repeated  here  itsing  the 
corner  diffraction  coefficients.  The  results  for  the  four  different  combinations  of 
incident  and  diffracted  polarizations  are  shown  in  Figures  14  to  17.  As  the  figures 
show,  the  previous  corner  diffraction  solution  is  discontinuous  at  the  false  shadow 
boundaries  {$  as  240°  and  $  as  300°)  of  the  front  and  back  edges. 


2.8  Stationary  Phase  Method 


The  stationary  phase  method  is  a  widely  used  method  for  asymptotically 
evaluating  a  certain  class  of  integrals  described  below.  The  results  are  stated. 
Further  informatior  on  the  stationary  phase  method  and  asymptotic  series  is  given 
in  [21]  for  the  following 

J{k)=  (2.17) 

«/a 


V 


Figure  13:  Defioition  of  angles  for  the  Corner  Diffraction  CoefHcien 
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Picvlout  Coiner  DifTraciiim 
Moment  Method 


Figure  14:  Co-polarized  RCS  in  the  <p  =  60°  plane  of  a  2A  square  plate  with  a 
polarized  fixed  source  at  6^  —  45°,  =  0°. 


Previoui  Corner  Diffraction 
Moroec!.  Method 


Figure  15:  Co-polsrizcd  RCS  in  the  4>  =  60°  plane  of  a  2A  square  plate  with  a 
polarized  fixed  source  at  3*  =  45°,  =  0°. 


0  60  120  eo  240  300  360 

0  (deg.) 

Figure  16:  Cross-polarized  RCS  in  the  <f>  —  60®  plane  of  a  2A  square  plate  with  a 
<9*  polarized  fixed  source  at  6^  =  45®,  (f>*  —  0®. 


Figure  17:  Cross-polarized  RCS  in  the  <f>  —  60®  plane  of  a  2A  square  plate  with  a 
polarized  fixed  source  at  9'  =  45®,  =  0®. 
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where  x  is  a  real  variable;  K  is  real,  positive  and  large;  ^(x)  is  a  real,  continuous 
function  with  continuous  derivatives  for  o  <  x  <  fe;  —  0,for  a  <  x$  <  b  and 

Xg  is  not  close  to  a  or  b.  Note  that  F(x)  is  a  complex  function  which  is  slowly 
varying  and  well-behaved  on  [a,b)  then  I{k)  is  asymptotically  approximated  by 


+  imeiNW-f)  -  imeiMW-51  +  O  (»-i)  .  (2.18) 

K(jy{a)  \  J 

The  first  term  is  the  contribution  from  the  stationary  point  while  the  second  and 
third  terms  are  the  contributions  from  the  end  points  b  and  a,  respectively.  If 
there  is  no  stationary  point  in  the  interval  [a,b]  and  there  are  not  any  stationary 
points  Xg  near  the  end  points  a  and  b,  the  integral  is  asymptotically  equal  to  the 
sum  of  the  end  point  such  that 

This  result  will  be  used  to  reduce  the  equivalent  edge  currents  to  comer  diffraction 
coefRcients.  It  is  also  easily  seen  that  the  terms  of  order  l//c  in  the  stationary 
phase  approximation  give  the  exact  integral  in  the  special  case  where  F(x)  is  a 
constant,  such  that  F(a)  =  F{b)  —  F.  Then,  ^(x)  =  cx  such  that  ^^(x)  =  c 
which  is  a  constant.  This  obviously  indicates  that  no  stationary  points  are  in  the 
interval  or  near  the  end  points.  Plugging  into  the  above  equation  and  retaining 
only  terms  of  order  I/k,  one  obtains  that 


/(k)  ~  ] 

K  C 


IF 

K  C 


(2.20) 


KJC  KJC 


(2.21) 
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v^hich  is  the  same  as  the  result  one  obtains  by  integrating  the  following  expression 

I{k)  -  .  {2.22) 
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CHAPTER  III 


THEORY 


3.1  Introduction 

This  chapter  outlines  the  derivation  of  the  new  equi'>alent  cuixents  for  finding 
the  far  zone  scattered  fields  from  a  perfectly  conducting  wedge  illuminated  by 
a  plane  wave.  This  is  a  combination  of  the  work  done  by  Michaeli  [13,14]  and 
Buyukdura  [22],  The  contribution  from  a  single  face  of  the  wedge  is  written  in  the 
form  of  an  equivalent  edge  current.  The  equivalent  currents  are  then  reduced  to 
corner  diffraction  coefficients.  Then  the  comer  diffraction  coefficients  are  written 
in  a  form  that  is  consistent  with  previous  forms  for  diffraction  coefficients. 

3.2  Physical  Optics  Surface  Integral 

Gordon  [23]  showed  that  the  far  zone  radiation  integral  over  the  PO  currents 
for  a  plane  wave  incident  on  a  flat  plate  may  be  converted  from  a  surface  integral  to 
a  line  integral  around  the  edge  of  the  scatterer.  Later  Buyukdura  [22]  arrived  at  the 
same  result  using  a  different  methoo.  The  steps  he  used  to  find  the  PO  equivalent 
edge  currents  are  repeated  het The  PO  equivalent  edge  current  is  foimd  for  a 
plate  with  the  0>fa;  s  illuminated.  The  results  may  be  easily  used  for  a  general 
wedge  since  the  FO  contribution  from  each  illuminated  face  may  be  calculated 
separately  and  s  tmmed  to  ^ve  the  total  PO  contribution.  If  a  face  of  the  wedge 
is  not  illumina'^ed  then  that  face  makes  no  contribution  to  the  scattered  field  and 
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does  not  interact  with  the  other  face  of  the  wedge  in  the  PO  approximation.  Using 
the  edge  fixed  coordinates  ^ven  in  Figure  1  and  choosing  the  normal  direction  so 
that  <  TT  it  is  assumed  that  the  contribution  from  the  0-face  to  the  far  zone 
scattered  field  may  be  written  as 

I  (3.1) 

where  c  is  the  curve  along  the  edge  or  edges  of  the  scatterer  and  X(<)  is  an  unknown 
function  of  the  source  and  observation  directions  which  is  to  be  found.  The  field 
produced  by  electric  and  magnetic  line  currents  placed  along  the  edge  are  given  by 

gs  =  f\ZQSXaxil  +  »xiM]  (3.2) 

An  R  Jc  ^  ■* 


where 

jkj.  =  ifei .  kz  —  ki'Z,  (3.3) 

Equating  this  result  with  Equation  (3.1)  one  finds  that 


I  —  Ki  -  JfjcotjScos^ 


(3.4) 


and 


where 


M  =  -ifjZo 


sin 

sin/d 


K{t)=iKi  +  bKi,, 


(3,5) 


(3.6) 


Placing  the  PO  currents  in  the  radiation  integral  over  the  surface  of  the  plate  gives 

^  X  j  X  /  /(2n  X  (3.7) 

4ir  it  J  Js 


where 

=  2n  X  H*(®,0,*)  (3.8) 
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(3.9) 

(3.10) 


fcj  =  ka'  •  X,  kg  =  ka^  •  z 

and  is  the  GO  incident  field.  Using  the  following  vector  identity 

J{kt-kl.)x  j{kz-kl)z 


ei(**-*i)*ei(**-*J)*  =  y .  V  X  I  [(ib;*  -  jb*)x  -  (ib*  -  jbi)z] 


jk^h^ 


with 


gives 


.2  (k.~kif+(k.-k\) 

= - p - 


i\2 


(3,11) 


(3.12) 


S’  =  ~  -  -J- i  X  i  X  /  (2ft  X  g{)!/  ■  V  X 

r  .  .  ,  ej(^x-ki)Xgj(kz-ki)z  1 

[{kz  -  kl)x  -  (fc*  -  fc*  )z] - JI^ - i  dxdz  .  (3.13) 


Then  applying  Stokes’  theorem,  one  obtains  that 


dz  (3.14) 


4ir  R 

where  the  line  integral  is  around  the  edge  of  the  plate.  Comparing  this  integral 
with  the  integral  in  Equation  (3.1)  results  in  the  following: 

}i _ if 

BJU  p  tUB  fp  T  Bin  f. 

H - 

jk 

!  -•  mn  /9  ro*  /A  4-  nin  /?’  rnn  th! 

(3.16) 


^  2  „i8in/3coB^  +  sin^'coB^^ 

Kt  =  - p - 

2  ^sin/Jcos^  +  sin/?^  cos^^ 

ff»  =  - p - 


(3.16) 


where 


using 


H'  =  hHi  +  nHi  +  iHi 


g*  =  i-ik*  X  # 
Zo 


(3.17) 


(3.18) 
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gives 


Hi  =  +  Hi  cot  /?'  cos  .  (3.19) 

ioQ  sin  p 

Using  Equations  (3.15),  (3.16),  and  (3.19)  in  Equations  (3.4)  and  (3.5)  gives 

M  =  ^^0  8  ^  U(ir  -  ^')  sin ^  cos  <f)  +  sin  0'  cos  <f>'  _ 

jk  *  sin  /?  ^(cos  —  cos /3^)^  +  (sin  jS  cos  <j)  +  sin  cos  )2 


and 


/  = 


X 

+ 


U(ir  -  ^') 

2Yo  sin  sin  /9  cos  ^  +  sin  /?'  cos 

jk  Sin/3'  (cos/S  —  cos  /3')^  +  (sin  /9  cos  ^  +  sin  j9'  cos  ^')^| 


2  i  (sin  p  cos  ^  +  sin/9'  cos  ^')(cot  /9  cos  ^  +  cot  /9'  cos  ^') 


jfc  *  (co8/9  -  cos/3')^  +  (sin  /9  cos  (j>  +  sin  /3'  cos  W 


'I 


(3.21) 


where  the  unit  step  function  U{v  —  ^')  has  been  included  to  stress  the  fact  that 
the  equivalent  currents  are  sero  if  the  face  is  not  illuminated  by  the  incident  Reid. 


Letting 


results  in 


sin  P  cos  <f>  ^  (cos  P  —  cos  P')^ 

Mn/3'  8in/9'(8in/9co8^  +  8in/3'cos^') 


^  2Zo  i  sin^  U{n-<p) 
jk  *  sin  P  sin  /?'  (cos  7  +  cos  4>') 


(3.22) 


(3.23) 


2Koj,^sin^'  Ujir-if,') 
jk  *  sin^/S'  (cos  7  +  cos  ) 

4.  2  ^jUj^  -4^)  (cot  jg cos ^  4-  cot P'  cos ^') 

jfc  *  sin^'  cos  7  4- cos 


(3.24) 


Since  the  above  equivalent  currents  (referred  to  later  as  the  LPO  component  of  the 
equivalent  currents  since  they  result  from  the  reduction  of  the  PO  surface  to  a  line 
integral)  were  derived  from  a  surface  integral  over  a  finite  surface  and  integrand, 
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it  might  be  expected  that  they  wotdd  remain  finite  for  all  aspects  of  incidence 
and  observation.  Even  though  the  currents  themselves  become  infinite  for  certain 
directions  of  incidence  and  observation,  the  fields  radiated  by  the  currents  remain 
finite,  as  explmned  below.  The  magnetic  current  becomes  infinite  as  sin  ^  — »  0. 
This  is  not  a  problem  however,  since  it  will  always  be  cancelled  by  a  factor  of  sin/3 
in  the  radiation  integral  due  to  the  i  x  t  term.  It  appears  that  the  electric  current 
becomes  infinite  for  edge  on  incidence  (sin/?'  =  0),  but  the  current  remains  finite 
here  since  the  factor  of  sin  /3'  in  the  denominator  is  cancelled  by  a  factor  of  sin  /?'  in 
El.  Both  of  the  currents  become  infinite  at  the  GTD  shadow  boundaries  (/8  =  /3' 
and  ^  ±  ^'  =  ±7r).  This  is  otdy  an  artifact  of  the  conversion  of  the  surface  integral 
to  the  line  integral.  When  the  complete  integration  is  done  avound  the  edges  of  a 
finite  scatterer,  the  infinite  contnbutions  from  the  different  edges  combine  to  yield 
a  finite  result.  This  is  illustrated  analytically  for  the  special  case  of  backscatter 
at  broadside  of  a  rectangular  flat  plate  in  Section  4.5  and  in  other  cases  by  many 
numerical  examples  in  Chapter  V. 

3.3  Fringe  Equivalent  Currents 

Many  people  have  worked  on  the  problem  of  finding  equivalent  edge  cur¬ 
rents  for  a  perfectly  conducting  wedge  that  are  valid  away  from  the  Keller  cone. 
Mitzner  [24]  derived  Incremental  Length  Diffraction  Coefficients  (ILDC)  using 
some  symmetry  arguments.  Knott  [25]  later  showed  that  Mitzner’s  ILDC  could 
be  written  as  eqtuvalent  edge  currents.  Michadi  [13]  later  derived  equivalent  edge 
currents  for  the  wedge  problem  using  a  more  mathematically  rigorous  method. 
Knott  [25]  also  showed  that  Mitzner’s  ILDC  and  Michaeli’s  equivalent  edge  cur¬ 
rents  are  the  same  except  that  Mitzner’s  ILDC  give  expressions  for  the  fringe  field 
only,  while  Michaeli’s  equivalent  edge  currents  are  for  the  total  field.  In  this  first 


37 


attempt  at  the  equivalent  currents,  the  plate  is  modeled  as  a  series  of  infinitesimal 
strips  oriented  normal  to  the  edge  (i.c.  the  strips  are  parallel  to  the  6  direction). 

In  a  later  paper  Michacli  [14]  points  out  many  singularities  in  the  original  ex¬ 
pressions  for  the  equivalent  fringe  currents  and  derives  new  expressions  which  have 
fewer  singularities.  This  is  accomplished  by  changing  the  direction  of  asymptotic 
evaluation  of  the  integrsd  from  normal  to  the  edge  to  a  direction  skewed  from  the 
edge  normal.  Not  only  does  the  new  orientation  of  the  strips  reduce  the  number  of 
singularities  in  the  fringe  equivalent  currents,  it  is,  as  will  be  explained  later,  also 
more  consistent  with  the  physics  of  the  problem.  The  original  expressions  for  the 
equivalent  fringe  currents  had  many  singrilarities.  Michacli  points  out  that  there 
are  so  many  singularities  that  it  is  probably  impossible  to  remove  all  of  them.  Go¬ 
ing  back  to  the  radiation  integral  over  one  face  of  the  wedge  he  points  out  the  cause 
of  two  of  the  singularities.  The  equivalent  fringe  currents  become  infinite  when 
s  -  S  =  5'  'S  where  S  is  the  observation  direction,  s'  is  the  incident  ray  direction,  and 
b  is  the  edge  fixed  coorcUnate  direction  defined  previously  in  Figure  1.  The  current 
becomes  infinite  here  because  the  phase  of  the  PO  component  of  the  surface  cur¬ 
rent  cancels  the  phase  of  the  outgoing  wave.  The  fringe  equivalent  currents  also 
become  infinite  when  5  •  S  ~  S  where  a  =  icos0'  -1-  Ssin/JV  The  fringe  currents 
are  infinite  here  because  the  phase  of  the  UTD  surface  current  cancels  the  phase  of 
the  outgoing  wave.  The  cone  of  rays  where  the  equivalent  currents  become  infinite 

A  * 

(s  •  6  =  d  •  &)  in  the  second  case  may  be  reduced  to  a  single  direction  by  a  change 
of  variables.  This  is  included  in  the  summary  below.  A  summary  of  the  procedure 
Michael!  [14]  used  to  derive  the  fringe  equivalent  edge  currents  is  given  below: 

1.  The  wedge  is  replaced  at  each  point  along  its  edge  by  an  infinite 
wedge  with  fiat  faces. 
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2.  The  currents  on  one  face  of  the  infinite  wedge  are  expressed  as  a  spectral 
integral  using  the  solution  to  the  two-dimensional  infinite  wedge  problem. 

3.  The  field  due  to  these  currents  is  expressed  as  a  triple  integral  (one  over  the 
spectrum  and  a  double  integral  over  the  surface  of  the  tangent  wedge), 

4.  The  order  of  integration  is  changed.  The  integration  over  the  spectrum  and 
the  integration  over  one  of  the  spatial  variables  (b)  are  interchanged. 

5.  The  integral  over  b  is  replaced  by  an  integration  along  the  a  direction  (see 
Figure  18),  and  the  integration  variable  b  is  replaced  by  <r  sin  in  the  inte¬ 
gral. 

6.  The  asymptotic  end  point  contribution  from  the  integral  over  the  first  spatial 
coordinate,  <7,  is  determined. 

7.  The  integration  over  the  spectrum  is  done,  after  the  path  has  been  distorted, 
using  the  method  of  steepest  descent.  The  contributions  from  the  poles  (PO 
contribution)  arc  removed  leaving  the  fringe  contribution  only. 

8.  The  remaining  integral  is  an  integration  along  the  edge  of  the  wedge. 

9.  The  radiation  integral  for  electric  and  magnetic  line  sources  along  the  edge 
is  compared  to  the  integral  remaining  in  8  and  the  equivalent  edge  currents 
are  identified  as: 


jf  _  jVTD  _  jPO 
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(3.25) 

(3.26) 


Figure  18:  Integr&Uoo  direction  (or  strip  orientation)  used  in  finding  the 

equivalent  ^nge  c*xrrents. 


where: 


= 

iVTD  ^ 


+ 


= 


I 


PO 


a 


= 


2jZs\n(i)  (l/n)8in[(7r  —  a)/n]c8ca 
k  sin  /3  sin  /3^  cos  [(tt  —  a)  /n]  —  cos  (<^V^)  ^ 


2i 


1/n 


Bin 


fc  sin CCS  {<!>' Jn)  —  cos  [{tt  -  a)  /n]  |  2  sin/?' 

sinfiV  —  a) /nl,  ^ 

- : - - — ^(cosacot/S  —  cot  ^  cos  > 

sin  a  J 

2icot/?'^j 

fcn  sin  /?'  * 


~2jZBm<l>U{7r-<i>') 


k  sin  /?  sin  ^'(cos  +  cos  a) 


2jU{ir  - 


pi 

Zsinfi'^^ 


k  sin  /3'(co8  +  cos  a) 

(cot  cos  <f>  4-  cot  ^  cos 

cos"^  fjt  =  —j  In  +  \//i2  - 
sin  ^  cos  (ff  (cos  0  —  cos  y9')  cos  /?' 


sin/?' 


sin^  /?' 


(3.27) 


(3.28) 

(3.29) 

(3.30) 

(3.31) 

(3.32) 


and  the  branch  cut  chosen  is  assumed  to  be  the  same  as  previously  obtained  by 
Michaeli  (13] 


-1  = 


J  j/l 


M<-1 
-1  <  M  <  1 

p>\ 


(3.33) 


Although  Michaeli  does  not  specify  the  branch  for  the  Ins  in  Equation  (3.31), 
it  seems  reasonable  to  take  the  principal  branch  of  the  natural  log.  Specifically 
In  z  =  In  |z|  +ja>rg(z)  where  0  <  arg{z)  <  ir.  From  the  equations  for  the  equivalent 
currents,  it  is  obvious  that  the  aTg{z)  may  also  be  taken  between  ir  and  2>r,  since 
the  equations  remiun  unchanged  if  a  is  replaced  by  2ir  —  a.  An  equivalent  method 
of  finding  a  is  to  choose  the  branch  of  cos'^  n  so  that  or  is  continuous  for  /u  =  ±1 
and  the  two  branches  are  chosen  to  run  from  a  =  Trtoa  =  ir  -  joo  and  a  =  0 
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Re  a 


H  fjL  increasing 


Figure  19:  Branch  for  a. 


to  a  =  0  —  joo  respectively,  as  shown  in  Figure  19.  Notice  that  the  last  term  on 
the  right  hand  side  of  the  expression  for  in  Equation  (3.28)  will  be  dropped 
since  it  will  be  cancelled  by  a  term  of  opposite  sign  due  to  the  other  face.  Since 
each  edge  is  made  up  of  two  faces  and  the  term  does  not  depend  on  which  face 
is  being  illuminated  this  term  will  always  be  cancelled  by  the  corresponding  term 
from  the  other  face  for  a  perfectly  conducting  wedge. 

The  only  difference  between  the  procedure  Michaeh  used  to  find  these  eqmva- 
lent  currents  and  the  procedure  he  used  to  find  his  previous  equivalent  currents  [13] 
is  in  step  5  listed  previously.  This  change  in  variables  not  only  results  in  fewer 
singularities  in  the  equivalent  currents,  but  also  makes  more  sense  physically  as 
Ufimtsev  points  out  in  (12).  The  field  on  one  of  the  plates  of  the  infinite  tangent 
wedge  should  only  depend  on  the  contribution  from  the  strip  on  which  the  field 
point  is  lying.  The  only  strip  orientation  whiclr  meets  this  condition  is  the  orienta¬ 
tion  described  step  5,  where  the  first  spatial  integration  is  taken  in  the  &  direction. 
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In  this  case,  all  points  along  the  strip  lie  on  the  Keller  cone  so  that  when  the  line 
integral  along  the  edge  is  evaluated  the  major  contribution  to  any  field  point  on 
the  strip  will  be  due  to  the  strip  itself.  All  other  directions  of  integration  result  in 
equivalent  currents  which  “contaminate”  the  field  along  their  strips  with  contri¬ 
butions  from  the  other  strips.  The  singularities  remaining  in  the  new  equivalent 
fringe  currents  arc  listed  below  : 

1.  The  currents  become  infinite  when  the  direction  of  observation  is  the  same 
as  the  direction  of  a  glancing  incident  ray,  that  is  the  forward  scattering 
direction  in  the  plane  of  the  plate  where  =  /3'  ,<i>  =  Q  ,  and  4''  — 
Michael!  refers  to  this  as  the  Ufuntsev  singularity. 

2.  and  do  not  tend  to  a  definite  limit,  but  remain  bounded  as  a  — »  d-. 

3.  Both  find  approach  finite  limits  for  edge  on  incidence  {fi'  — ♦  0  or 
P'  tt),  except  of  course  when  s  =  which  corresponds  to  the  Ufimtsev 
singularity  described  in  item  1. 

4.  Both  and  are  discontinuous  for  4’  = 

5.  and  Af-^  are  finite  for  i*d’  =  or  ai:4’  =  ’Tj  although  the  components 
of  the  currents  become  infinite  here. 

Michaeli  points  out  that  all  of  these  singularities  except  for  the  first  one  are  in- 
iegrable  so  the  fringe  fields  foimd  using  these  currents  are  finite  for  all  aspects 
of  incidence  and  observation  except  in  a  single  forward  scatter  direction.  Since 
the  Ufimtsev  siugularity  is  in  the  forward  scattering  direction,  it  will  probably  not 
cause  any  trouble  in  solutions  of  practical  importance.  These  singularities  may 
cause  spikes  and  discontinuities  in  the  pattern  when  the  equivalent  currents  are 
used  as  the  bases  for  the  new  eorner  di&actioo  coefficient.  This  behavim  will  be 
illustrated  in  Ch^>ter  IV. 
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3.4  Equivalent  Currents  Written  in  the  Form  of  Diffraction  Coeffi¬ 
cients 

The  equivalent  edge  currents  may  be  written  in  a  diffraction  coefficient  form 
closer  to  that  used  by  Kouyoumjian  and  Pathak  [3].  Using  the  following  trigono¬ 
metric  identities: 


1 

cos  a  -f-  cos  b 


1 

cos  a  -f  cos  h 


4- 


sin  a 

cos  a  —  cos  b 
sin  6 

cos  a  -  cos  6 


and  a  little  algebra,  one  obtuns  the  equivalent  edge  currents  (0-face  contribution 
only)  given  by 


j^LPO  ^  j^VTD  _  j^PO 

(3.38) 

le  + 

(3.39) 

rLPO  ,  rVTD  rPO 
■*e,m  *e,m 

(3.40) 

where: 


= 

Ck(7),C2h.#') 

jLPO 

•  < 

2 

c. 

.  J 
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MUTDjUTD 

jUTD 


1  C|i(a),C2(a,ir-a) 


•  '  M^POjLPo' 

[PO 

jLPO 

LPO  jLPO 

”*  (with  7  replaced  by  o)  (3.41) 


ChiS)  = 


.  c  _  \  /3  421 

Kq^  sin/? sin/?' sin  5  ’  *  ZQk  sin^/?' 


^ \  cotjScoB^  +  cot|9'cos£ 

=  ■* - - 


COS  7  = 


cos  a  = 


jb  8injsin;9' 

sin/? cos ^  (cos^-coB/?')^  ..v 

sin^'  Bin/?'(8in/?coB^  + sin^'cos^') 
sin  /?  cos  ^  (cos  /?  -  cos  /?')  cos  /?' 
sin/?'  sin^^' 


(3.43) 


(3.45) 

(3.46) 


•  0  ,ir-^'<0 

U*  =  (3,46) 

1  ,ir  -  >  0 

The  branch  of  coe~^  is  described  in  Section  3.3.  For  the  special  case  of  a  flat  plate 
(n=2),  it  is  much  easier  to  consider  the  contribution  from  each  edge  due  to  both 
faces  than  to  consider  the  contribution  from  the  edge  due  to  each  face  separately. 
For  the  case  of  a  flat  plate  the  total  contribution  from  each  edge  is  given  by  (note: 
the  normal  may  be  taken  in  either  direction  perpendicular  to  the  plane  of  the 
plate.): 

U  =  -  M '’O  (3.47) 

/  =  /.  +  /m  (3.48) 

Ajn  =  +  {3.«) 
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where: 


jLPO 
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(p«) 


±  cot 

-v. 

1 

+ 

1 

1 

lUTD 

2 

c. 

I  |eot  +  cot 

|cot  ^!Lz(|±iO  j  +  „i  ^,r±(°,+  ^')  j'  I 


MLPOjIPO  ■ 

jlPO 

5*  = 


(3.50) 


{with  7  replaced  by  a) 


-1  ,ir-^'<0 

1  yn-(f/>0 

where  C2,  C*,  7,  and  a  are  the  same  as  the  expressions  given  earlier  in  Equa¬ 
tions  (3.42)  to  (3.45).  The  terms  may  be  associated  with  the  0-face  incident 
and  reflection  shadow  boundaries  and  the  N-face  incident  and  reflection  shadow 
boundaries,  as  indicated  previously  for  the  UTD  diffraction  coefficients. 

3.6  New  Corner  Diffraction  Coefficient 


Buyukdura  (22)  has  derived  the  corner  diffraction  coefficients  from  the  equiv¬ 
alent  currents  using  the  stationary  phase  approxinoation.  If  the  object  is  made 
up  of  flat  faces  with  straight  edges  then  ail  of  the  edges  are  terminated  in  sharp 
vertices.  It  is  assumed  that  the  current  near  eadi  edge,  even  very  near  the  vertex, 
is  unperturbed  by  the  other  edge  making  up  the  vertex.  The  same  assumption 
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was  made  in  the  derivation  of  the  previous  corner  diffraction  coefficient  given  in 
Section  2.7.  Despite  this  shortcoming  the  corner  diffraction  solution  gives  fairly 
accurate  results  (see  the  next  chapter  for  numerical  results).  Since  it  is  assumed 
that  the  two  edges  forming  the  vertex  do  not  interact  near  the  vertex  (to  first 
order)  the  scattered  field  may  be  expressed  as  a  line  integral  along  the  edges  of  the 
object.  The  line  integral  along  the  edges  may  obviously  be  expressed  as  the  sum 
of  integrals  along  each  individual  edge  such  that 

E^  =  y\y, [Zo  i  X  i  X  t  /  +  i  X  t  M]  (3.51) 

^  Air  R  Join  ^ 

where  the  first  sum  is  over  the  faces  of  the  object,  the  second  sum  is  over  the  edges 
of  each  face,  and  the  edge  fixed  coordinates  I  smd  b  along  with  the  limits  on  the 
integrals  obviously  depend  on  which  face  and  edge  are  being  considered.  Looking 
only  at  the  contril>ation  due  to  one  face  and  one  edge  at  a  time,  one  finds  that 


E% n)  =  [Zo  «  X  i  X  <  J-  +  a  X  i  m]  (3.52) 

where  the  coordinates  are  defined  in  Figure  20.  Notice  that  along  the  straight  edge 
the  edge  fixed  coordinate  directions  (^,n,t),  and  the  angular  directions  ()9',^', 
are  constant.  The  only  factors  in  the  equivalent  currents  which  vary  along  the 
edge,  therefore,  are  and  Rf.  To  simplify  the  integration,  these  factors  will  be 
separated  out  using  the  following  definitions: 


kninfinnP’  ksin/^sin/?' 


M  = 


2jl  ZoGm^t 


(3.53) 

(3.54) 


fe  sin  ^  sin 

From  which  it  is  easily  seen  that  Gim,  and  Gm  are  all  constant  along  the 
straight  edge.  Taking  the  phase  reference  at  some  arbitrary  point  0  along  the 
t-axis  and  using: 


Ei(z)  = 


(3.56) 


47 


(3.56) 


a;  =  0 

(3.57) 

kz  =  k  cos 

(3.58) 

gives 

S’-  -  '  ‘T  /I* X 3 X i  f !(“>  + 

27r  R  [sm/Ssm/S'  sin/JsmjS' 

sm/jsin/j'  j 

(3.59) 

Since  the  term  in  brackets  is  constant,  using  either  the  method  of  stationary  phase 

or  regular  integration  (sec  Chapter  II)  and  only  including  the  contribution  from 

the  end  point  at  b  along  with  the  following  relations: 

fit  ”  ^0 

(3.60) 

(3.61) 

i  xi  ~  — ^siuj9 

(3.62) 

s  X  i  X  <  =  fisinP 

(3.83) 

results  in  the  following  expressions  for  the  new  diffraction  coefficients  in  terms  of 
the  new  equivalent  currents: 
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(3.64) 


Di  D^2 

^-jha 

0  Dl 

.4. 

a 

►  • 

=  ± 

+ 

1  sin  0  sin  0^ 
in  cos  0  —  cos  0‘ 

hlYoE* 

MIZqHI 

.^2. 

Im/Hi 

(3.66) 


where  and  £7^  are  the  0  and  ^  components  of  the  corner  diffraction  contribution 
from  one  corner,  one  edge,  and  one  face.  The  phase  is  referenced  to  the  comer. 

and  JS^i  are  the  0'  and  components  of  the  incident  field  at  the  corner. 
The  angles  are  the  edge  fixed  angles  with  the  origin  located  at  the  corner  as 
illustrated  in  Figure  21.  The  equivalent  currents  /e>  -fm*  and  M  are  given  by 
Equations  (3.41)  and  (3.39),  respectively,  or  by  Equations  (3.49)  and  (3.47)  for 
the  case  of  a  flat  plate.  The  sign  used  depends  on  which  end  point  the  particular 
contribution  is  coming  from  and  the  direction  of  integration.  The  direction  of 
integration  is  always  chosen  to  be  in  the  counterclockwise  direction  when  the  face 
is  viewed  from  the  positive  n-direction.  The  negative  sign  is  used  for  the  first 
corner  contribution  of  each  corner  (as  the  direction  of  integration  is  traced  out) 
and  the  positive  sign  is  used  for  the  second  comer  contribution  of  the  edge  as 
indicated  in  Figure  22.  Equations  for  the  comer  diffraction  coefficients  with  the 
expressions  given  in  Equations  (3.41)  and  (3.39)  substituted  into  Equation  (3.65) 
are  given  in  Appendix  A. 


S.6  Conclusion 


Expressions  for  new  equivalent  currents  based  on  the  PTD  are  given  in  Equa¬ 
tions  (3.39)  and  (3.41).  The  new  equivalent  currents  are  used  to  derive  a  new 
comer  diffraction  coefficient.  Both  the  new  equivalent  currents  and  the  new  cor- 
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Edge  of 
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Figure  21:  Definition  of  the  Angles  used  in  the  New  Corner  Diffraction 

Coefficients. 
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Direction  of  Integration 


^  V  coordinates  .  , 

XI  shown  for  the  front  oottom  edge) 

Figure  22:  Direction  of  Integration  for  one  face. 

ner  diffraction  pve  the  same  resxilts  for  the  first  order  scattering  from  a  convex 
body  made  up  of  perfectly  conducting  flat  plates.  Even  though  the  equivalent  cur¬ 
rents  are  not  valid  in  the  region  of  a  vertex  (corner),  it  will  be  shown  in  Chapter  V 
that  the  equivalent  currents  (and  the  corner  diffraction  coefficients)  give  reason¬ 
able  results  for  many  directions  of  incidence  and  observation.  Several  important 
properties  of  tlus  solution  are  discussed  in  Chapter  IV  followed  by  several  examples 
in  Chapter  V. 
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CHAPTER  IV 


PROPERTIES  OF  THE  NEW  CORNER  DIFFRACTION 

SOLUTION 


4.1  Introduction 

The  behavior  of  the  new  corner  diffraction  coefficient  will  be  discussed  for 
various  situations.  The  significance  of  its  various  properties  will  be  outlined  to 
^ve  insight  into  how  the  solution  may  be  expected  to  behave  in  certun  special 
cases.  Several  more  examples  given  in  Chapter  V  illustrate  that  the  solution  agrees 
closely  with  either  the  Method  of  Moments  or  previous  first  order  solutions  in  most 
cases. 

4.2  Singularities 

First  the  singularities  in  the  solution  will  be  discussed.  The  singularities  in  the 
new  corner  diffraction  coefficients  are  simply  the  combination  of  the  singularities  in 
the  line  PO  component  and  in  the  firinge  component  of  the  new  equivalent  currents. 
The  singularities  and  some  examples  illustrating  how  they  may  effect  a  pattern  are 
given  below. 

1.  Ufimtsev’s  singularity  in  the  forward  scattering  direction  in  the  half  plane 
of  the  face  (/3  •—  ,^  =  0  ,  and  4"'  =  ’>’)•  This  singularity  was  previously 

discussed  in  item  1  of  Section  3.3.  To  illustraie  how  this  may  affect  the 
results,  the  bistatic  RCS  of  a  two  wavelength  square  plate  is  cslculated.  In  the 
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Figure  23:  PUte  Geometry  for  the  Dluatration  of  Uhirntrev’s  Singularity. 

firut  example  the  source  remains  hxed  at  9^  =  90°  and  =  0°  (see  Figure  23) 
while  the  pattern  is  taken  in  the  ^  =  0°  plane.  The  results  for  the  H-plane 
pattern  in  the  x-e  plane  are  shown  compared  with  a  Method  of  Moments 
stdution  in  Figure  24.  This  result  shows  that  the  new  comer  diffraction 
solution  is  fairly  accurate  for  directions  greater  than  15°  to  20°  away  from 
the  forward  direction  for  cases  of  grasing  incidence.  The  previous  corner 
diffraction  solution  and  equivalent  current  solution  also  become  infinite  in  the 
forward  direction  for  grasing  inddence.  The  previous  solutions  are  compared 
with  the  new  solution  in  Figure  25  for  the  !T-plane  pattern  and  the  geometry 
^ven  above.  The  new  and  old  solutions  are  actually  so  close  that  they 
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0  &C  120  180  240  300  360 

0  (deg.) 

Figure  24:  H-plane  pattern  for  <f>  =  0°  cut  of  a  2A  square  plate  with  a  fixed 

source  at  6i  =  90®,^*  =  0®. 


Figure  25:  H-plane  pattern  for  <f>  =  0°  cut  of  a  2A  square  plate  with  a  fixed 

source  at  6x  =  90°,  =  0°. 
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0  60  120  180  240  300  360 

G  (deg.) 

Figure  26:  H-plaae  pattern  for  (f>  =  45°  cut  of  a  2A  square  plate  with  a  fixed 

source  at  =  90°,  =  45°. 

are  indistinguishable  in  the  plot  in  Figure  25.  To  see  how  moving  off  the 
Keller  cone  affects  the  behavior  of  the  solution,  the  bistatic  RCS  of  the  two 
wavelength  plate  is  plotted  for  a  pattern  in  the  ^  -  45°  plane  wliile  the  source 
is  held  fixed  at  =  90°  and  =  45°.  The  results  for  the  H-plane  pattern 
are  plotted  along  with  Moment  Method  calculations  in  Figure  26.  In  this 
example  the  new  comer  diffraction  solution  is  fairly  accurate  in  the  major 
lobe  region  away  from  the  forward  direction  {9  —  90°).  The  previous  corner 
diffraction  solution  and  previous  equivalent  current  solution  are  compared 
with  the  new  solution  for  this  example  in  Figure  27.  Although  the  pattern 
is  taken  away  from  the  Keller  cone,  all  three  solutions  are  fairly  close. 

2.  iL  (and  D2)  and  (and  do  not  tend  to  definite  limits  as  s  -♦  0^ ,  where 
7  —  tsin^  +  5cos/3'  and  is  shown  in  Figure  18,  but  they  remain  bounded. 
In  practice  this  means  that  both  D2  and  17^,  and  therefore  and  are 


0  60  120  ttO  240  300  360 

0  (deg.) 


Figure  27:  H-plane  pattern  for  ^  =  45°  cut  of  a  2 A  square  plate  with  a  fixed 

source  at  Bi  —  90°, =  45°. 


discontinuous  at  this  point  in  the  pattern.  A  simple  example  illustrates  how 
this  discontinuity  can  affect  a  pattern.  The  bistatic  RCS  from  a  flat  plate, 
two  wavelengths  on  a  side,  in  the  x-y  plane  is  considered.  The  source,  linearly 
polarised  in  the  direction,  remains  flxed  at  9*  —  45°  and  =  0°  while 
the  pattern  is  taken  near  the  x-y  plane  (5  =  89°).  The  source  geometry 
is  shown  in  Figure  28.  The  bistatic  RCS  is  given  in  Figures  29  and  30  for 
the  co-polarised  and  cross  polarized  fields,  respectively.  The  abrupt  null  at 
fti  135°  in  the  co-polarized  pattern  and  the  spike  at  the  same  location  in 
the  cross-polarised  pattern  are  due  to  discontinuities  in  the  contribution  from 
edge  4  (indicated  in  Figure  28).  The  point  ^  s::  135°  coincides  with  ~ 
and  ^4  ft;  0  where  ^4,  ^4,  and  ^4  are  the  edge  fixed  coordinates  for  edge  4. 
Due  to  the  geometry  ^4  ft;  ^  so  the  discontinuity  in  is  due  to 

the  discontinuity  in  Z>2  and,  likewise,  the  discontinuity  in  <r^  is  due  to  the 
discontinuity  in  D^. 
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Figure  28:  2A  square  plale  in  the  x*y  plane  illuminated  by  a  fixed  somce  located 

in  the  6*  —  45®,  =  O'*  direction. 
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I  I  1  I  I  I  —I  I  I  ,  1  I 

0  30  60  so  120  eo  «0 

♦  (deg.) 

Figure  29:  RCS  for  the  6  =  89®  cut  of  a  2A  square  plate  with  a  polarized  fixed 

source  at  6*  —  45®,  =  0®. 


Figure  30:  RCS  for  the  9  -  89®  cut  of  a  2A  square  plate  with  a  polarized  fixed 

source  at  9*  =  45®,  =  0®. 
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Figure  31:  Fi^,  0'  =  136® 

The  specific  terms  which  axe  discontinuous  in  and  are  examined  fur¬ 
ther  here.  The  discontinuity  in  at  this  point  is  due  to  the  discontinuity 
in  Specifically,  it  is  due  to  the  following  factor: 

cos  a  cot  0*  —  cot  0  cos  ^ 


FUPJ.'t')  = 


Sin  a 


(4.1) 


■^im  is  plotted  as  a  function  of  0  with  ^  as  a  parameter  and  0'  =  135®  in 
Figure  31.  For  ^  on  the  order  of  a  few  degrees,  Fim  changes  sign  abruptly  at 
0  —  0\  This  is  exactly  what  happens  to  the  contribution  to  the  $  polarised 
field  from  edge  4  in  the  above  example.  The  same  conclusion  may  be  reached 
analytically  by  reahsing  that 

cos  a  cot  0'  -  col  0  cos  ^  ^0 


FMP’A) 

where 


sin  a 


- - — (4.2) 

^(A^8inj9')2  +  (A/?)2 

^  (4.3) 
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=  )8  -  /?'  (4.4) 

for  A<{>  and  Aj9  small.  Which  will  further  reduce  to  the  following  signum 
function  if  A(f>  «  A/?: 

P  ta  at  ±\  COS  a  cot /3' -  cot cos  .... 

P  ,  <f>)  - - - - - -  sgn{Afi)  (4.5 

sin  a 

as  <j>  -*  0. 


The  discontinuity  in  is  examined  next.  The  only  factor  in  that  is  dis¬ 
continuous  at  the  intersection  of  the  Keller  cone  and  the  half  plane  associated 
with  the  edge  (/3  =  and  ^  =  0)  is 


P  _  Bin  4> 

"m  —  • 

sma 

This  factor  only  causes  a  discontinuity  in  the  term  since  another 

factor  in  the  causes  to  go  to  zero  here.  Fm  is  plotted  as  a 

function  of  0  with  ^  as  a  parameter  and  ~  135°  in  Figure  32.  For  ^ 
on  the  order  of  a  few  degrees,  Fm  becomes  a  very  narrow  spike  at  =13', 
This  also  occurs  to  the  contribution  for  the  6  polarized  field  from  edge  4  in 
the  above  example.  The  same  conclusion  may  be  reached  analytically  by 
realizing  that: 


where: 


A<t>  =  4> 


(4.7) 


(4.8) 


A3 = ^  -  /j* 


(4.9) 


for  A^  and  A/9  small.  Equations  (4.2)  and  (4.7)  show  that  Fjm  <^d  Fm  ere 
very  similar  with  the  major  difference  being  that  the  roles  of  A/9  and  A^ 
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Figure  32:  Fm,  /3'  =  135'^ 

axe  interchanged.  Fim  and  Fm  are  plotted  again  in  Figure  33  and  Figure  34, 
reapectively,  this  time  as  a  function  of  ^  with  ^  as  a  parameter  and  0'  =  135^. 


These  discontinuities  are  not  restricted  to  the  case  of  bistatic  scattering  and 
may  also  appear  for  backscatter  examples.  The  backscattered  field  (co-  and 
cross-polarised),  for  a  B  pdarised  field  incident  upon  the  plate  shown  in  Fig¬ 
ure  28  and  calculated  close  to  the  x-y  plane  {6  =  89*’),  is  shown  in  Figure  35. 
In  both  cases  the  patterns  have  a  glitch,  due  to  the  contributions  from  the  far 
edge  at  ^  =  90°,  similar  to  those  in  the  previous  bistatic  results.  It  should 
be  noted  that  the  fields  for  a  ^  polarised  incident  field  (shown  in  Figure  36) 
are  smooth  for  this  pattern.  The  cross  polarised  field,  <7^,  is  not  visible  on 
this  {dot  since  it  is  below  -50  through  the  entire  pattern. 
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Co-polviied  RCS 
Crow’polariied  RCS 


The  discontinuity  in  the  new  equivalent  currents  (and  the  diffraction  coef¬ 
ficients)  at  the  intersection  of  the  Keller  cone  and  the  infinite  half  plane 
associated  with  the  edge  and  ^  =  0)  may  be  expected  to  cause 

discontinuities  or  narrow  spikes  depending  on  the  polarization  and  the  pat¬ 
tern  cut.  As  the  examples  illustrate  these  disturbances  only  affect  a  typical 
pattern  cut  for  around  5®  to  10®.  Further  information  on  how  the  solution 
behaves  in  this  region  is  discussed  in  Section  4.3. 

3.  Excluding  the  case  when  the  diffracted  ray  is  in  the  same  direction  as  the 
incident  ray  (see  item  1),  the  currents  are  all  finite  for  both  edge  on  incidence 
{fi'  —*  0,ir)  and  edge  on  observation  (/?  — »  0,7r).  This  is  shown  for  edge 
on  incidence  in  Appendix  B.  For  edge  on  observation  it  is  easily  seen  by 
noting  that  the  factors  of  l/sin/3  in  the  equivalent  currents  are  cancelled  by 
a  factor  of  sin  0  in  the  numerator  of  the  radiation  integral  due  to  either  sxi 
or  s  X  3  X  t  depending  on  whether  it  is  an  equivalent  electric  current  or  an 
equivalent  magnetic  current. 

4.  D2  wid  are  discontinuous  for  =  tt,  since 

are  discontinuous  for  <f>'  =  w,  except  for  two  important  special  cases.  In  the 
case  of  backscatter  ^  and  all  of  the  components  with  the  factor 

contiuning  the  step  function  are  zero  for  =  tt,  so  that  the  resulting  corner 
diffraction  coefficient  is  continuous.  When  the  observation  point  is  on  the 
Keller  cone  and  so  that  the  discontinuous 

components  in  the  corner  diffraction  coefficients  cancel  leaving  the  corner 
diffraction  coeffidents  continuous.  It  is  pointed  out  how  this  discontinuity 
may  affect  some  typical  patterns  in  Section  4.4. 
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6.  Dg,  Z?2>  *“^<1  infinite  at  the  shadow  boundaries  because  the  LPO 

components  are  infinite.  The  LPO  components  are  infinite  because  the  vec¬ 
tor  identity  used  to  derive  them  is  not  valid  at  the  shadow  boundaries  (see 
the  end  of  Section  3.2).  This  is  not  a  problem  in  practice  since  the  cor¬ 
ner  diffraction  coefficients  will  remain  finite  a  small  distance  away  from  the 
shadow  boundaries. 

6.  All  of  the  diffraction  coefficients  become  infinite  on  the  Keller  cone.  This 
is  not  a  problem  however,  once  for  straight  edges  the  contribution  from  the 
other  comer  of  the  edge  will  also  be  infinite  and  their  sum  will  combine  to 
form  the  familiar  on  x  jx  pattern  associated  with  a  constant  current  element. 

7.  The  contribution  to  the  scattered  field  from  the  individual  corner  happens 
to  be  infinite  when  0  =  0'.  This  is  not  the  case  for  the  finite  plate,  however, 
since  the  individual  contributions  from  all  of  the  edges  balance  each  other 
to  give  a  bounded  result.  This  is  shown  analytically  in  Section  4.5  for  the 
special  case  of  backscatter  from  a  rectangular  plat  at  normal  incidence.  It 
is  also  illustrated  many  times  in  Chapter  V  since  0  =  0'  (or  at  least  one  edge 
in  ail  of  the  patterns  taken  in  the  principle  plane. 

4.S  Boundary  Conditions 

Since  the  new  solution  is  only  valid  in  the  far  zone  it  does  not  need  to  meet 
the  boundary  conditions  on  the  scattering  structure,  except  for  the  special  case 
of  a  half  plane  (semi-infinite).  Actually  the  scattered  field  in  the  plane  of  a  given 
face  must  meet  different  conditions.  The  conditions  on  the  scattered  field  may  be 
easily  found  using  the  far  zone  radiation  integral  combined  with  the  fact  that  the 
faces  of  the  structure  are  flat.  Given  an  arbitrary  perfectly  conducting  flat  face 
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Observation 

Direction 


shown  in  Figure  37  and  illuminated  from  some  arbitrary  direction,  the  far  zone 
scattered  field  is  given  by 

jkZQe-i^^. 


4ir  R 


xsx  I  izdy  (4.10) 


where  the  important  coordinates  are  defined  in  Figure  37  and  J#  is  the  actual 
current  on  the  surface  of  the  face.  Obviously  Jg  only  has  nonzero  components  in 
the  X  and  y  directions.  As  a  result  when  the  observation  point  is  in  the  plane 
of  the  face  (the  x-y  plane)  the  z  component  of  the  total  scattered  electric  field 
must  be  zero.  This  also  implies,  by  reciprocity,  that  the  total  scattered  field  in 
the  special  case  when  the  face  shown  is  the  only  scattering  object  (n=2)  should  be 
aero  everywhere  for  a  i  directed  source  placed  in  the  x-y  pla.ne.  These  conditions 
only  apply  to  the  total  scattered  field  and  do  not  necessarily  apply  to  a  first  order 
solution  as  b  given  here. 

It  is  assumed  that  the  first  order  fields  radiated  by  the  equivalent  currents 
must  meet  these  conditions  in  the  plane  of  the  face.  (Why  this  assumption  is  not 
correct  will  be  discussed  later  in  this  section.)  This  assumption  combined  with  the 
radiation  integrd  for  the  equivalent  edge  currents  given  by 

1 1^0  a  X  i  X  <  /  +  ^  X  <  Af]  (4.11) 

or  the  comer  diffraction  coefficients  lead  to  the  following  properties: 


1.  M=0  and  =  0  when  =  0  or  ^  =  ir  , 


2.  /m  =  Af  =:  0  and  DJ  =  ^  when  —  0  and  n  =  2  or  =  ir  and  n  =  2 

(/«  =:  0  and  Dg  contributes  zero  to  the  scattered  field  since  =  0) 

Pl*\gging  ^  =  0,^  into  the  new  solution  it  is  seen  that  the  first  property  is  met 
everywhere  in  the  plane  of  the  face  with  the  exception  of  two  directions.  One 
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is  in  the  forward  scatter  direction  which  corresponds  to  the  Ufimtsev  singuiaritv 
discussed  in  item  1  of  Section  4.2,  The  other  direction  where  this  property  is  not 
met  is  the  direction  where  the  Keller  cone  and  the  infinite  half  plane  associated 
with  the  edge  meet.  It  is  interesting  to  note  that  this  direction  lies  alon^  the 
infinite  strip  used  in  deriving  the  equivalent  fringe  currents.  T\.o  of  the  equivalent 
currents  associated  with  the  solution  to  the  infinite  wedge  aad  ^  are 

discontinuous  here,  as  illustrated  in  item  2  of  Section  4.2.  It  appears  that  the 
solution  meets  the  first  far  zone  condition  given  above  everywhere  except  on  this 
strip  where  it  meets  the  boundary  conditions  on  a  PEC.  Plugging  in  =  0,7r  and 
n  =  2  into  Equations  (3.47)  and  (3.49)  it  is  easily  seen  that  the  second  condition 
given  above  is  not  met  in  general. 

Both  the  previous  equivalent  currents  given  in  Section  2.6  and  the  previous 
corner  diffraction  solution  given  in  Section  2.7  meet  the  boundary  conditions  on 
the  infinite  half  plane  associated  with  the  edge.  That  is  /  =  0  (and  Dl  =  0)  for 
^  =  0  and  7  =  0  (and  Dg  =  0)  for  (f/  =  0. 

It  has  been  shown  that,  in  general,  the  first  order  solution  should  not  meet 
the  far  zone  conditions  described  here.  Specilically,  both  Sikta  [26]  (for  the  special 
case  of  backscatter)  and  Tiberio  [27]  have  shown  that  for  the  two-dimcnsioneii 
problem  of  an  infinite  strip  of  finite  width  the  double  diffraction  term  cancels 
the  first  o’-der  term  so  that  the  far  zone  conditions  in  the  plane  of  the  face  are 
met.  \'he  new  solution  is  still  consistent  with  this  requirement  since  it  reduces 
to  the  previous  equivalent  currents  for  observation  points  on  the  Keller  cone  (see 
Section  4.6).  Since  the  solution  to  the  infinite  quarter  plane  has  not  been  reduced 
to  a  ray  optical  form,  the  conditions  that  the  first  order  field  must  meet  in  the 
plane  of  the  face  when  the  observation  point  is  not  on  the  Keller  cone  arc  not  clear. 
In  this  case  the  new  solution  and  the  previous  solutions  behave  much  differently. 
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Therefore,  when  the  observation  point  is  not  on  the  Keller  cone  the  higher  order 
terms  for  the  new  solution  may  be  significantly  different  than  those  developed  for 
previous  solutions.  Some  results  for  three-dimensional  problems  including  higher 
order  terms  have  been  published  by  Sikta,  et  al  (2|  and  [26]  and  Michaeli  [28]  using 
the  previous  equivalent  currents  and  the  new  equivalent  currents,  respectively. 

4.4  Reciprocity 

Although  only  the  total  scattered  field  must  satisfy  reciprocity,  once  again  it 
would  seem  desirable  for  the  first  order  scattered  field  given  by  the  new  equivalent 
currents  or  corner  diffraction  coefficients  to  also  satisfy  reciprocity.  Both  the  re¬ 
flected  GO  field  and  the  edge  diffracted  field  satisfy  reciprocity  even  though  they 
only  make  up  a  partial  contribution  to  the  total  scattered  field.  Even  though  the 
new  solution  is  not  reciprocal,  several  examples  are  given  that  show  it  is  probably 
close  enough,  in  most  cases,  for  engineering  purposes. 

A  simple  example  is  used  to  illustrate  that  the  new  solution  is  not  reciprocal. 
This  is  shown  by  considering  where  the  singularities  of  the  fringe  occur.  Neither 
Im  —  nor  approach  definite  limits  as  the 

observation  point  approaches  the  intersection  of  the  Keller  cone  and  the  infinite 
half  plane  associated  with  the  edge  of  the  plate.  However  if  the  source  is  allowed  to 
approach  this  position  in  the  same  way  all  of  the  currents  are  finite  and  approach 
definite  limits.  Using  these  properties  a  simple  example  is  given  which  shows  that 
the  solution  is  not  reciprocal.  A  triangular  flat  plate,  as  shown  in  Figure  38,  is 
illuminated  by  a  source  positioned  above  the  plane  of  the  plate  in  position  one. 
The  observation  point  is  placed  at  position  two  in  the  plane  of  the  plate  such  that 
it  lies  on  the  Keller  cone  for  the  points  on  edge  one,  as  indicated  in  Figure  38. 
Both  ill  and  are  undefined  for  points  along  edge  one  which  makes  a  nonzero 
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Position  1 


V 


contribution  to  the  radiation  integ;rai.  Therefore  the  scattered  field  is  not  defined 
here.  However  if  the  source  and  receiver  are  interchanged  all  of  the  currents  are 
well  defined  and  therefore  the  scattered  field  will  be  defined.  It  is  not  surprising 
that  the  new  solution  is  not  reciprocal  nnce  it  contains  PO  terms  which  are  not 
reciprocal  themselves. 

Several  examples  of  bistatic  scattering  from  a  flat  plate  or  a  cube  are  given 
here  to  illustrate  that  the  solution  is  essentially  reciprocal  in  many  instances.  All 
of  the  examples  axe  in  terms  of  linearly  polarised  antennas  in  the  far  zone  of  the 
scatterer.  First  the  RCS  is  calculated  when  the  first  antenna  is  receiving  and  the 
second  antenna  is  transmitting.  The  RCS  is  then  calculated  for  the  reciprocal 
problem  (antenna  one  transmitting  and  antenna  two  receiving).  Both  patterns  are 
plotted  on  the  same  graph  as  a  function  of  the  position  of  antenna  one,  so  that 
they  may  be  easily  compared. 

The  first  set  of  examples  show  the  scattering  from  a  square  plate  two  wave¬ 
lengths  on  a  side  as  shown  Figure  39.  In  the  first  two  examples  both  antennas  are 
moved  in  the  principle  plane  while  they  are  offset  by  90®  in  that  plane,  as  shown  in 
Figure  39.  Figure  40  shows  the  results  when  both  antennas  are  linearly  polarized 

A 

in  the  9  direction,  while  Figure  41  shows  the  results  when  they  are  both  linea.dy 
polarized  in  the  ^  direction.  The  cross-polarized  fields  are  zero  in  both  cases.  The 
results  show  that  in  this  case  the  solution  is  reciprocal  for  practical  purposes.  This 
is  not  surprising  since  in  this  case  the  major  contribution  to  the  fields  comes  from 
the  front  and  back  edges.  The  receiving  antenna  lies  on  the  Keller  cone  of  the 
front  and  back  edges  for  the  entire  pattern  cut.  The  new  solution  is  reciprocal 
for  observation  points  on  the  Keller  cone  since,  as  is  shown  later  in  Section  4.6,  it 
reduces  to  Keller’s  result  there. 

In  the  next  two  examples  both  antennas  are  moved  in  the  ^  =  45®  plane  offset 
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Figure  40:  RCS  of  2A  square  plate  for  9  polarized  antennas  (which  are  positioned 
as  indicated  in  Figure  39)  in  the  ^  =  0'^  plane. 
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Figure  42:  Relative  positions  of  the  two  antennas  for  the  next  two  patterns. 

by  90°  in  that  plane,  iw  shown  in  Figure  42.  Figure  43  shows  the  results  when 
both  antennas  are  linearly  polarized  in  the  $  direction,  while  Figure  44  shows  the 
results  when  they  are  both  linearly  polarized  in  the  ^  direction.  The  cross-polarized 
fields  are  zero  in  both  cases.  Although  the  results  are  not  nearly  as  close  as  in  the 
previous  patterns,  the  major  differences  occur  when  either  the  receiver  or  the  source 
is  near  the  plane  of  the  plate  where  higher  order  terms  may  become  important. 
The  antennas  are  moved  in  the  ^  =  30°  plane  offset  by  30°  in  that  plane  as  shown 
in  Figure  45  for  the  next  three  examples.  Figure  46  shows  the  results  when  both 
antennas  are  linearly  polarized  in  the  9  direction,  while  Figure  47  shows  the  results 
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Figure  43;  RCS  of  2A  square  plate  for  9  polarised  antermaa  (which  are  positioned 
aa  indicated  m  Figure  42)  in  the  <j>  ~  45°  plane. 


Figure  44:  RCS  of  2A  square  plate  for  i  polarised  antennas  (which  are  positioned 
as  indicated  m  Figure  42)  in  the  ^  =  45°  plane. 
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Figure  46:  RCS  of  2A  square  plate  for  9  polarized  antennas  (which  are  positioned 
as  indicated  in  Figure  45)  in  the  ^  =  30°  plane. 

when  they  are  both  linearly  polarized  in  the  ^  direction.  In  this  case,  however, 
the  cross-polarized  fields  are  no  longer  zero.  It  should  be  noted  that  reciprocity 
does  not  imply  that  =  <r^  for  bistatic  scattering,  and  in  general  ^ 
for  bistatic  scattering.  Reciprodty  implies  that  the  pattern  remains  unchanged  if 
the  antenna  originally  transmitting  receives  while  the  antenna  origin^ly  receiving 
transmits,  but  the  positions  of  the  two  antennas  remain  unchanged  (in  the  bistatic 
case,  give  the  results  when  the  positions  of  the  antennas  are  interchanged). 

Figure  48  gives  the  cross-polarized  results  when  antenna  one  is  linearly  polarized 
in  the  9  direction  and  antenna  two  is  linearly  polarized  in  the  ^  direction.  If 
the  polarizations  are  reversed,  that  is  antenna  one  is  linearly  polarized  in  the  4> 
directive  and  antenna  two  is  polarized  in  the  9  direction  one  simply  obtains  the 
mirror  image  of  both  curves  due  to  the  symmetry  of  the  problem. 

The  co-polaiized  fields  are  close  to  being  reciprocal  in  this  case.  The  cross- 
polarized  fields  are  reasonably  reciprocal  for  repons  away  from  the  plane  of  the 
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Figure  48:  RCS  of  2A  square  plate  for  cross  polarized  antennas  (which  are 
positioned  as  in^cated  in  Figure  45)  in  the  ^  =  30”  plane. 
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Figure  49:  Relative  positions  of  the  two  antennas  for  the  next  four  patterns. 

plate.  In  regions  near  the  plane  of  the  plate,  the  pattern  with  the  null  at  ^  =  90° 
is  probably  more  accurate  since  the  plate  does  not  radiate  any  6  pciarized  field  in 
this  plane.  The  lack  of  a  null  here  in  the  other  case  may  be  explained  by  the  fact 
that  Im  (or  D2)  does  not  meet  the  conditions  in  the  plane  of  the  face  described  in 
the  previous  section.  Finally  the  antennas  are  moved  in  a  plane  found  by  rotating 
the  x-s  plane  120°  in  the  ^  direction  followed  by  tilting  it  30°  in  the  d  direction, 
as  shown  in  Figure  49.  The  two  antennas  are  then  offset  by  150°  in  the  x-y  plane. 

Figure  50  shows  the  results  when  both  antennas  are  linearly  polarised  in  the 
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direction  parallel  to  the  plane  the  pattern  is  taken  in  and  perpendicular  to  the 
incident  ray  (Horizontal  Polarization),  while  Figure  51  shows  the  results  when  they 
are  both  linearly  polarized  in  the  direction  perpendicular  to  the  plane  the  pattern 
is  taken  in  (Vertical  Polarization).  Once  again  the  cross-polarized  fields  are  not 
zero.  Figure  52  gives  the  cross-polarized  results  when  antenna  one  is  horizontally 
polarized  and  antenna  two  is  vertically  polarized.  Figure  53  gives  the  patterns 
for  the  opposite  case  where  antenna  one  is  vertically  polarized  and  antenna  two  is 
horizontally  polarized.  The  discontinuity  in  the  four  patterns  which  occurs  when 
4>  =  120®  and  antenna  two  is  transmitting  (solid  curves  in  Figures  50,  51,  52, 
and  53)  is  due  to  the  discontinuity  in  the  PO  components  described  in  item  4  of 
Section  4.2.  The  discontinuity  in  the  four  patterns  which  occurs  when  4>  =  90®  and 
antenna  two  is  transmitting  is  due  to  the  discontinuity  in  the  contribution  from 
the  edge  along  the  x-axis.  The  direction  ~  90®  corresponds  to  the  intersection  of 
the  Keller  cone  and  the  infinite  half  plane  associated  with  this  edge  {0i  %  %  30® 


81 


90  t20 

4>  (deg.) 


Figure  52:  RCS  of  2A  square  plate  for  antenna  one  horisonially  polarized  and 
antenna  two  vertically  pmariced  (which  are  positioned  as  indicated  in  Figure  49 


which  are  positioned  as  indicated  in  Figure  49). 
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Figure  53:  RCS  of  2A  square  plate  for  antenna  one  vertically  polariaed  and 
antenna  two  horizontaliy  polarized  (which  are  positioned  as  indicated  in 

Figure  49). 

and  <{>i  w  0°)  where  the  comer  diffraction  coefficient  is  discoutinuous  (see  item  2 
of  Section  4.2).  The  discontinuity  in  the  four  patterns  which  occurs  when  ^  =  90“ 
and  antenna  one  is  trsmsmitting  (dashed  line  in  Figures  50,  51,  52,  and  53)  is  due 
to  the  discontinuity  in  the  PO  components  described  in  item  4  of  Section  4.2.  For 
directions  at  least  15“  to  20“  away  from  these  discontinuities  the  solution  is  dose 
to  being  redprocal. 

To  get  some  indication  of  the  effect  of  wedge  angle  on  these  examples  the 
same  patterns  are  repeated  for  the  case  of  a  cube  two  wavelengths  on  a  side.  The 
cube  is  positioned  in  the  usual  way  so  that  each  edge  is  parallel  to  one  of  the  three 
axes.  The  relative  podtions  of  the  antennsa  for  the  first  two  patterns  ere  given  in 
Figure  39.  Figure  54  shows  the  results  when  both  antennas  arc  linearly  polarized 
in  the  8  direction,  while  Figure  56  shows  the  results  when  they  arc  both  linearly 
polarized  in  the  ^  direction.  As  in  the  case  of  the  fiat  plate  the  co'polarized  field 
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Figure  55:  RCS  of  2A  cube  for  6  polarized  antennas  ^hich  are  positioned  as 
indicated  in  Figure  42)  in  the  ^  =  45°  plane. 

in  the  principle  plane  is  reciprocal  for  practical  purposes.  This  is  expected  since 
the  major  contribution  to  the  scattered  field  is  from  the  edges  perpendicular  to  the 
pattern  plane.  In  both  examples  the  receiving  antenna  lies  on  the  Keller  cone  of 
these  edges  where  the  solution  reduces  to  Keller’s  form  as  discussed  in  Section  4.6. 
The  small  discontinxiity  at  the  peaks  of  the  ^  polarized  pattern  is  due  to  the 
contribution  from  the  two  edges  parallel  to  the  pattern  plane  and  farthest  from 
the  source  and  receiver.  The  contribution  from  these  two  edges  is  discontinuous 
in  this  case  due  to  the  discontinuity  in  and  described  earlier  in  item  2  of 
Section  4.2. 

In  the  next  two  examples  both  antennas  are  moved  in  the  4>  =  45°  plane  offset 

by  90°  in  that  plane,  as  shown  in  Figure  42.  Figure  56  shows  the  results  when  both 

« 

antennas  are  linearly  polarized  in  the  0  direction,  while  Figure  57  shows  the  results 
when  they  are  both  linearly  polarized  in  the  ^  direction.  The  cross-polarized  fields 
are  zero  in  both  cases.  The  results  are  furly  close  to  being  reciprocal,  although 
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Figure  57:  RCS  of  2A  cube  for  ^  polarined  antenn&s  (which  are  positioned  as 
indicated  in  Figure  42)  in  the  ^  =  45^  plane. 

not  as  close  as  in  the  principle  plane.  All  of  the  patterns  are  discontinuous  at 
—  0®, 90®,  180®.  This  discontinuity  is  due  to  the  discontinuity  in  the  LPO  and 
PO  terms  of  and  as  described  in  item  4  of  Section  4.2.  The  antennas  are 
moved  in  the  ^  =  30®  plane  offset  by  30°  in  that  plane  as  shown  in  Figure  45  for  the 
next  three  examples.  Figure  58  shows  the  results  when  both  antennas  are  linearly 
polarised  in  the  6  direction,  while  Figure  59  shows  the  results  when  they  are  both 
linearly  polarised  in  the  ^  direction.  In  this  case  the  cross-polarised  fields  are  not 
sero.  In  order  to  avoid  numerical  problenos  associated  with  edge  on  incidence  (see 
Appendix  B)  the  cube  was  tilted  5®  in  the  plane  normal  to  the  ^  =  30®  plane  for 
these  calculations.  Figure  48  gives  the  croes-polarised  results  when  antenna  one 
is  linearly  pdarised  in  the  $  direction  and  antenna  two  is  linearly  polarised  in  the 
^  direction.  If  the  polarisations  are  reversed,  that  is,  if  antenna  one  is  linearly 
polarised  in  the  ^  direction  and  antenna  two  is  polarised  in  the  $  direction,  one 
simply  obtains  the  mirror  image  of  both  curves  due  to  the  symmetry  of  the  problem. 
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The  co-poiarised  fields  are  essentially  reciprocal  in  this  case.  The  cross-polarized 
fields  are  close  to  being  reciprocal,  except  in  the  region  from  6  =  60°  to  90°.  All 
of  the  patterns  for  the  case  corresponding  to  the  final  example  given  for  the  plate 
have  so  many  discontinuities  that  it  is  very  difficult  to  get  any  information  from 
them.  For  this  reason  they  have  been  omitted  here. 

The  major  difference  between  the  results  for  the  flat  plate  and  the  results  for 
the  cube  is  the  increase  in  the  number  of  regions  in  space  where  the  PO  disconti¬ 
nuity  (item  4  of  Section  4.2)  shows  up  in  the  patterns. 

The  examples  given  here  illustrate  that  the  solution  is  close  ‘  u  1 nf  tt^ciprocal 
for  the  cases  shown,  except  in  regions  of  space  where  the  patterns  are  discontinuous. 

4.6  RCS  of  a  Rectangular  Plate  at  Broadside 

As  a  check  on  the  corner  diffraction  coefficient,  the  radar  cross  section  of  a 
flat  rectangular  plate  (length  i,  width  v;,  and  area  A)  at  broadside  was  cuiculated 
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analytically  (see  Appendix  C).  The  resulting  co-polarized  cross  section,  for  either 
polarization,  is  given  by 

This  is  the  well  known  PO  result,  which  is  generally  accepted  as  adequate  for  engi¬ 
neering  purposes  for  sufficiently  large  plates.  Therefore  the  new  corner  diffraction 
solution  can  be  expected  to  give  the  correct  value  for  backscatter  from  a  general 
shape  (consisting  of  flat  plates)  at  broadside. 

4.6  Equivalent  Currents  on  the  Keller  Cone 


The  expressions  for  the  equivalent  currents  on  the  Keller  cone  are  easily  ob- 
tsdned  by  letting  =  /3Mn  Equations  (3.38)  and  (3.40)  combined  with  the  fact 
that  a  =  7  =  ^  in  this  case.  The  resulting  expressions  for  the  equivalent  currents 
with  the  contribution  from  both  faces  included  is: 


M  =  - 


+ 


/  =  - 


(4.14) 


These  are  the  same  expressions  for  the  equivalent  currents  as  obtained  by  Ryan 
and  Peters  (1]  given  in  Section  2.6.  Ryan  and  Peters  equivalent  currents  predict 
the  same  fields  as  Keller  [10]  obtained  if  the  edges  are  straight  and  the  observation 
point  is  on  the  Keller  cone.  Since  the  new  equivalent  currents  and  the  new  corner 
diffraction  coefficient  give  the  same  result,  the  corner  diffraction  coefficient  reduces 
to  Keller’s  result  for  observation  points  on  the  Keller  cone. 
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4.7  Physical  Significance  of  Parameters  a  and  7 


The  two  angles  a  and  7  in  the  expressions  for  the  eqtiivalent  currents  obvi¬ 
ously  do  not  represent  physical  angles  since  they  may  become  complex  for  several 
directions  of  incidence  and  observation.  It  is  easily  seen,  however,  that  they  are 
related  to  the  directions  sJong  which  the  asymptotic  endpoint  contribution  is  taken 
in  deriving  the  equivalent  currents  (described  for  the  equivalent  fringe  currents  in 
item  5  of  Section  3.3).  The  first  term,  in  Equations  (3.44)  and  (3.45)  is  common 
to  both  a  and  7  along  with  the  previous  a  parameter  obtained  by  Michaeli  [13]  in 
his  first  paper  on  equivalent  currents.  The  second  term  depends  on  the  direction 
of  integration  used  in  finding  the  equivalent  currents.  Buyukdura  [29]  has  shown 
that  the  equivalent  currents  associated  with  the  PO  surface  integral  (LPO  terms) 
may  also  be  obtained  by  making  two  simple  changes  to  the  procedure  Michaeli 
used  to  find  the  fringe  equivalent  currents.  The  GO  current  is  used  on  the  infinite 
tangent  wedge  instead  of  the  true  current  and  the  spatial  integral  is  taken  in  the 
direction  such  that 


cot  5 1  =  "T 


cos  (3  —  cos 


(4.15) 


sin^flcos^-b  wny3'  cos^' 

V  ”  ^  ;8  the  angle  between  the  positive  edge  tangent  and  the  direction  of  in¬ 
tegration  as  shown  in  Figure  61.  As  explained  in  Chapter  III  the  spatial  integral 
used  to  evaluate  the  fringe  equivalent  currents  was  taken  in  a  direction  62 1  where 
cot92  =  cot/?',  also  shown  in  Figure  61. 

Using  the  above  angles  it  has  been  found  that  the  new  equivalent  currents  (and 
corner  diffraction  coefficients)  may  be  written  in  a  symmetric  form  for  the  special 
case  of  a  flat  plate.  In  tlus  symmetric  form,  the  arguments  of  the  trigonometric 
functions  are  real  valued  for  all  angles  of  incidence  and  observation.  It  is  also 
speculated  that  this  may  be  done  for  a  general  wedge  angle,  although  a  method  of 
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functions  are  real  valued  for  all  angles  of  incidence  and  observation.  It  is  also 
speculated  that  this  may  be  done  for  a  general  wedge  angle,  although  a  method  of 
doing  so  has  not  been  found.  In  the  special  case  of  a  flat  plate  (n=2)  the  currents 
may  be  rewritten  by  properly  using  Equations  (3.47)  and  (3.49)  along  with  the 
trigonometric  identities  given  in  Section  3.4.  The  resulting  expressions  for  the  new 
equivalent  currents  are  given  by 


Af  =  ^  mUTD  _  mPO 


I  =  le  +  l 


m 


T  -  r^PO  ,  rUTD  tPO 
■*e,m  —  ie,m  •+*  ^e,m  ~  -‘e, 


m 


(4.16) 

(4.17) 

(4.18) 


where: 
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C/l  = 


Cj  — 


C2  =  - 


a2(m)  =  ^  ^  ^ 

k  sm^m  ^  ' 

jHi  sin^  /3'  sin  ^  cos  ^ 

^'b^  sin  0  sin  V>2  sin  ^ 

sin  sin  ^ 

^0^  sin  0'  sin  ^2  cos  ^ 

C2  =  -  ^  ^  |cot  )9'  (cos  V>2  -  cos^  i8') 

*  sin  ^2  sin 

-  sin^/S*  cot/Scos^l 

co8V>m  =  fiin$msin/3cos^  +  cos^mcos^ 

cosV^ln  =  sinflmsin/?^  cos^’ -  coB®rnC08/3* 

_  COS0~COS0^ 

sin/Scos^  +  sin/3*  cos^' 
cot  ^2  =  cot  /3^ 

S‘  =  I  •’^-^'<0 

I  1  ,ir  -  >  0 


(4.20) 


(4.21) 


(4.22) 


S'  = 


(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 


(4.28) 


The  details  are  outUned  in  Appendix  D.  Notice  that  the  arguments  of  the  cotan¬ 
gent  functions  now  correspond  to  physical  angles  and  will  always  be  real  valued. 
The  an^es  V*li  i^2^  ^2  ^  angles  between  either  the  incident  or  the 

diffracted  ray  and  either  of  the  spatial  integration  directions,  6}  or  62,  used  in 
finding  the  equivalent  cunents.  The  four  new  angles  are  defined  by: 


cos^l  =  i-d^i 


cos^ 


j  •  ©2 


.>  • 
s  ’Oy 


(4.29) 

(4.30) 

(4.31) 


cosV»2 
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where: 


d'l  =  Ssin^i -f  f  CO8  01 
a’2  =  b  sin  ^2  +  ^  cos  $2 


(4.33) 

(4.34) 


and  and  62  are  given  in  Equations  (4.26)  and  (4.27),  respectively.  The  first 
two  angles,  rpi  and  are  the  angles  between  the  observation  direction  and  the 
direction  that  the  asymptotic  endpnipt  contribution  was  taken  for  the  LPO  and 
fringe  components,  respectively.  Likewise  the  second  two  angles,  xj)'^  and  ^2)  are 
the  angles  between  the  incident  ray  direction  and  the  direction  that  the  asymptotic 
endpoint  contribution  was  taken  for  the  LPO  and  fringe  components,  respectively. 
The  LPO  and  PO  components  are  easily  written  in  this  way.  The  UTD  components 
look  more  complicated  in  this  form,  however  it  may  be  possible  to  write  the  UTD 
components  in  a  form  similar  to  the  LPO  and  PO  components  if  the  endpoint 
contribution  is  done  over  entirely  in  a  different  coordinate  system. 
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CHAPTER  V 


RESULTS 


5.1  Introduction 

The  scattered  fields  predicted  by  the  new  corner  diffraction  coefficient  derived 
in  Chapter  III  are  compared  with  previous  high  frequency  solutions,  measurements, 
and  Moment  Method  solutions  for  some  specific  examples.  The  new  solution  is  first 
compared  with  a  few  backscatter  examples  commonly  found  in  the  literature  to 
show  that  in  the  principal  planes  it  generally  gives  the  same  results  as  previous 
solutions.  Then  the  new  corner  diffraction  coefficient  is  compared  with  measure¬ 
ments  for  a  swept  frequency  case  which  is  described  later.  Finally  the  new  corner 
diffraction  solution  is  compared  with  Method  of  Moments  and  measurements  for 
some  new  examples.  Throughout  this  chapter  the  new  corner  diffraction  solution 
and  the  new  equivalent  current  solution  will  both  be  referred  to  as  the  new  or  new 
corner  solution.  Note  that  it  was  shown  in  Chapter  III  that  both  give  the  same 
reauKc  for  first  order  scattering. 

5  .2  Brief  Description  of  Program 

A  computer  program  written  previously  at  OSU  to  calculate  the  RCS  of  a 
convex  perfectly  conducting  body  constructed  from  flat  plates  with  strmght  edges 
was  modified  to  calculate  the  RCS  using  the  new  solution.  Since  the  program 
was  written  to  handle  objects  with  straight  edges,  the  equivalent  current  form 
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of  the  solution,  given  in  Equations  (3.38)  and  (3.40),  was  used.  The  integration 
along  the  edge  was  done  in  closed  form  to  give  the  contribution  from  each  edge. 
For  structures  other  than  flat  plates,  it  is  also  necessary  to  take  shadowing  into 
account.  Since  only  concave  structures  are  considered  here  this  is  fadrly  simple. 
The  shadowing  of  the  incident  field  is  accounted  for  in  the  and  terms  in  the 
equivalent  currents  and  the  unit  step  functions  in  the  LPO  and  PO  components  of 
the  equivalent  currents.  The  shadowing  of  the  diffracted  ray  is  more  complicated. 
The  contribution  from  the  LPO  and  PO  components  are  present  everywhere.  They 
are  treated  as  if  they  are  radiating  in  free  space,  just  like  typical  Physical  Optics 
surface  currents.  The  UTD  components  are  shadowed  like  diffracted  fields.  They 
do  not  contribute  if  the  observation  point  is  inside  the  wedge  >  nir). 

5.3  Examples 

5.3.1  Ross 

The  backscatter  at  9.227  GKz  from  a  4”  x  4”  plate  in  the  y-z  plane  is  shown 
next  (sec  Figure  62)  with  the  results  for  a  pattern  cut  taken  in  the  x*z  plane  are 
shown  in  Figure  63,  Both  the  new  corner  and  Sikta’s  results  are  basically  the  same 
near  broadside,  as  expected.  Sikta’s  results  include  double  and  triple  diffraction, 
and  thus  are  more  accurate.  In  Figure  64  the  new  solution  is  compared  with  the 
equivalent  currents  of  Ryan  and  Peters  (the  first  order  scdution  used  by  Sikta)  and 
the  previous  corner  solution.  The  three  first  order  methods  give  essentially  the 
same  answer  in  this  case. 

6.3.2  Northrop  Fin 

The  backscatter  from  the  Northrop  fin  shape  shown  in  Figure  65  is  compared 
to  measurements  and  calculations  made  by  Sikta  (26).  The  E- plane  pattern  in 
the  x-y  plane  is  shown  at  17.76  GHz  in  Figure  66.  The  new  solution  is  compared 
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Figure  63:  E-pla^ie  pattern  in  X'Z  plane  of  Robs  plate 
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Figure  65:  Nortlux>p  Fin. 
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Figure  66;  E-plane  pattern  in  x-y  plane  at  17.76  GHz  of  Northrop  Fin. 

with  measurements  [30],  Sikta’s  results,  and  results  obtained  by  Northrop  [19]. 
Sikta’s  results  include  higher  order  terms  and  are  closer  to  the  measured  results 
than  the  new  corner  diffraction  solution  which  includes  only  first  order  terms.  The 
results  obtained  by  Northrop  [19],  which  are  based  on  the  exact  solution  of  an 
infinite  strip,  are  also  close  to  the  measured  results  and  are  essentially  the  same 
as  Si*  .r/s  results.  The  new  solution  once  again  agrees  well  with  the  measurements 
in  th  j  re^on  near  broadside  to  the  plate.  The  new  solution  is  compared  to  first 
order  results  obtained  by  Sikta  in  Figure  67,  to  show  that  higher  order  terms  are 
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important  in  regions  farther  than  30®  from  broadside.  Once  agedn  Sikta  used  the 
Ryan  and  Peters  equivalent  currents  to  obtain  the  first  order  contribution  to  the 
scattered  field.  The  new  solution  and  the  first  order  solution  obtained  by  Sikta  are 
about  the  same  as  expected. 

6.3.3  Sikta’s  triangle 

The  new  solution  is  first  compared  with  calculations  and  measurements  ob- 
tedned  by  Sikta  [26].  The  backscatter  in  the  principal  plane  from  the  triangle  shown 
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Figure  68:  Sikta’s  Triangle. 


in  Figure  68  was  calculated  for  different  values  of  a/A  and  a.  The  fields  calculated 
using  the  new  corner  diffraction  coefficients  are  compared  with  the  results  given  by 
Sikta  [26].  Sikta  calculated  the  fields  using  Ryan  and  Peters  equivalent  currents  to 
find  the  first  order  terms  plus  higher  order  terms  (double  and  triple  diffraction). 
The  results  are  shown  in  Figures  69  to  92.  Throughout  these  examples  the  two 
different  polarizations  are  alternated,  with  the  6  polarized  field  being  followed  by 
the  ^  polarized  field.  In  Figures  69  to  T4,  the  triangle  length  is  varied  (o  =  2A, 
3 A,  9A)  while  the  vertex  angle,  a,  is  held  constant  at  30°. 


Figure  09:  The  E-plane  pattern  {<f>  =  0°)  for  Sikta’s  triangle  with  o  =  2A  and 

a  =  30®. 
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Figure  70:  The  H-plane  pattern  (0  =  0®)  for  Sikta’s  triangle  with  o  =  2A  and 

a  =  30®. 


In  Figures  75  to  80  the  results  obtained  using  Ryan  and  Peters  [1]  equiva¬ 
lent  edge  currents  (the  first  order  solution  used  by  Sikta)  are  compared  with  the 
results  obtained  using  the  new  corner  diffraction  coefficient.  The  different  polar¬ 
izations  and  triangle  sizes  are  given  in  the  same  order  as  that  used  previously  in 
Figures  69  to  74.  In  Figures  81  to  88  the  triangle  length  is  held  constant  at  a  ==  4A 
while  the  vertex  angle  a  is  varied  (a  =  45®,  60®,  90®). 
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Figure  VI:  The  E-pl»ne  pattern  =  0®)  for  Siktn’e  triangle  with  a  =  3A  and 

a  =  30®, 
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Figure  72:  The  K'plane  petlern  {4>  —  0*’)  for  Sikta’s  triangle  with  a  =  3A  and 

a  =  30^ 
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Figure  74:  The  H-pUne  pattern  (^  =  0°)  for  Sikta’s  triangle  with  a  =  9A  and 

Q  =  30®. 
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Figure  75:  The  Ei-plane  pattern  =  0®)  for  Sikta’s  triangle  with  a  =  2A  and 

a  =  30®. 
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Figure  76:  The  H'plane  pattern  =  0°)  for  Sikta’s  triangle  with  a  =  2A  and 

a  =  30°. 
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Figure  77;  The  E-pI&ne  pattern 
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=  0°)  for  Sikta’s  triangle  with  a  =  3A  and 
o  =  30”. 
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Figure  78:  The  H-plune  pattern  (^  =  0®)  for  Sikta’s  triangle  with  o  =  3A  and 

a  =  30®. 
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Figure  79:  The  E-plane  pattern  {<f>  =  0®)  for  Sikta’s  triangle  with  a  =  9A  and 

a  =  30®. 
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Figure  81:  The  E-plane  pattern  =  0®)  for  Sikta’s  triangle  with  a  =  4A  and 

a  =  45®. 
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Figure  82:  The  H-pIane  p&ttern 


120 


60 


90 

0  (deg.) 


I 

CO 


=  0°)  for  Sikta's  triangle  with  a  =  4A  and 
a  =  45®. 
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Figure  83:  The  E-plane  pattern  {(f>  =  0®)  for  Sikta’s  triangle  with  a  =  4A  and 

a  =  60®. 
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Figure  84:  The  H-plane  pattern  =  0°)  fc 
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Figure  85:  The  E^pl&ne  p&ttera  (^  —  O'*)  for  Sikta’e  triangle  with  a  =  4A  and 

o  =  90^ 
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Figure  86:  The  H-plaoe  pattern  =  0®)  for  Sikia’s  triangle  with  a  =  4A  and 

a  =  90°. 
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before,  the  results  obtained  using  Ryan  and  Peters  [1]  equivalent  edge 
currents  (the  first  order  solution  used  by  Sikta)  are  compared  with  the  results 
obtained  using  the  new  corner  diffraction  coefficient  in  Figures  87  to  92.  The 
different  polarizations  and  triangle  sizes  are  given  in  the  same  order  as  that  used 
previously  in  Figures  81  to  86. 


Figure  87:  The  &pl&ne  pattern  {4>  =  0“)  for  Sikta's 

a  =  45“. 
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Figure  88:  The  H-plane 


i  0®)  for  Sikta's  triangle  with  o  =  4A  and 
a  p=  45®. 
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Figure  89:  The  E-pl&ne  pattern  =  0*^)  for  Sikta’s  triangle  with  a  =  4A  and 

a  =  60®. 
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Figure  90:  The  H*plane  pattern  (^  =  0“)  for  Sikta’s  triangle  with  o  -  4A 

a  =  60®. 
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Figure  91:  The  Emplane  pattern  =  0®)  for  Sikta’s  triangle  with  o  =  4A  and 

a  =  90®. 
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In  ^  cases  the  new  sdution  and  the  measurements  match  well  in  the  region 
within  50°  of  the  broadmde  to  the  plate.  In  addition  for  the  case  of  ^  polarization 
the  new  solution  and  the  measurements  are  close  enough  for  engineering  purposes 
all  the  way  to  the  plane  of  the  plate  in  the  region  away  from  the  tip  of  the  triangle. 
For  either  polarization  the  new  comer  and  the  measurements  differ  significantly 
in  the  re^on  of  the  tip  of  the  triangle.  Sikta  shows  that  higher  order  terms  are 
important  in  this  region.  The  new  solution  and  the  first  order  terms  obtained  by 
Sikta  using  the  equivalent  currents  of  Ryan  and  Peters  are  approximately  the  same 
for  all  of  the  examples  given  here.  There  are  only  minor  differences  between  the 
two  different  methods  which  become  smaller  as  the  size  of  the  triangle  is  increased. 

6.3.4  Cube 

A  few  examples  are  illustrated  here  for  the  scattering  from  a  cube  with  one 
wavelength  long  sides.  The  results  are  compared  with  those  obtained  using  the 
Method  of  Moments.  The  cube  is  positioned  relative  to  the  axes  in  the  usual  way 
(each  edge  is  parallel  to  one  of  the  coordinate  axes).  The  results  for  co-polarized 
backscatter  in  the  principle  plane  are  shown  for  the  H-plane  (with  9  =  90°)  in 
Figure  93  and  for  the  Emplane  (with  ^  =  90°)  in  Figure  94.  The  results  obtained 
are  very  close  to  those  obtained  using  the  Method  of  Moments. 

5.3.5  6”  Cube 

In  this  example  the  new  solution  is  compared  to  backscatter  measurements 
made  on  a  6”  cube.  All  of  the  results  are  for  co-polarized  fields  only.  With 
the  cube  oriented  so  that  each  edge  is  parallel  to  one  of  the  axes  as  shown  in 
Figure  95,  the  pattern  was  taken  in  the  x-y  plane  at  10  GHz.  The  results  for  a 
vertically  polarized  (or  9  polarized)  incident  field  are  shown  in  Figure  96,  while 
the  results  for  a  horizontally  polarized  (^  polarized)  incident  field  are  shown  in 


130 


Corner  Diffraction 
Moment  Method 
Normalised  to  10.7  dBfX* 


Figure  93:  H-plane  pattern  for  Cube  (with  9  =  90®). 
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Figure  94;  E- plane  pattern  for  Cube  (with  =  90°  )• 
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Figure  95:  Cube  geometry  with  the  illuminated  corners  labeled. 

Figure  97.  The  results  agree  fairly  well  everywhere,  with  differences  of  only  a  few 
dB  throughout  the  patterns.  The  patterns  are  repeated  for  a  cube  tilted  45°  in 
the  x-z  plane  as  illustrated  in  Figure  98.  The  results  for  the  H>pianc  and  E-plane 
patterns  taken  in  the  x-y  plane  are  given  in  Figures  99  and  100,  respectively.  The 
results  do  not  agree  as  well  in  this  case.  This  is  probably  due  to  a  combination  of 
errors  in  the  theory  and  in  the  measurements.  The  major  error  in  the  theory  is  the 
lack  of  higher  order  terms.  The  error  in  the  measurements  is  likely  two  fold.  First 
the  faces  of  the  cube  were  misaligned  and  so  they  did  not  form  edges  as  sharp  as 
may  be  required.  Secondly  it  seems  that  there  was  some  deviation  from  the  desired 
pattern  cuts  as  can  be  seen  from  the  lack  of  symmetry  in  the  measured  patterns. 

The  final  patterns  are  for  a  cube  tilted  45°  in  the  x-z  plane  followed  by  tilting 
it  35.2°  in  the  y-z  plane  so  that  the  z-axis  coincides  with  one  of  the  long  diagonals 
passing  through  the  center  of  the  cube  (as  illustrated  in  Figure  101).  The  results 
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for  the  H-plane  and  E-plane  patterns  taken  in  the  x-y  plane  are  given  in  Figures  102 
and  103,  respectively.  As  with  the  previous  patterns  the  higher  level  lobes  agree 
fairly  well,  but  the  lower  level  lobes  do  not  agree  well.  The  reasons  for  this  are  the 
same  as  those  stated  previously. 

5.S.6  Swept  Frequency  Results 

In  this  section  swept  frequency  measurements  ano  calculations  are  trans¬ 
formed  to  the  time  domain  as  described  by  Dominek  [31]  and  Leeper  (32].  A 
brief  description  of  the  idea  behind  this  is  given  here.  While  the  target  remains 
stationary  the  RCS  is  measured  over  a  wide  frequency  range.  An  Inverse  Fourier 
Transform  is  then  preformed  on  the  frequency  domain  data  to  give  the  ‘^finite 
bandwidth  impulse  response”  of  the  target  as  a  function  of  time.  In  this  manner 
the  different  scattering  centers  on  a  target  may  be  resolved  due  to  their  different 
electrical  lengths  from  the  source.  Obviously  the  closer  the  scattering  center  is 
to  the  source  the  earlier  its  contribution  will  appear  in  the  time  response.  The 
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Figure  103:  E-piane  p&twerc  for  a  6”  Cube  tilted  45°  in  the  x-8  plane, 

in  the  y-z  plane. 
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resolution  of  the  scattering  centers  in  the  time  domain  is  approximately  equal  to 
the  reciprocal  of  the  bandwidth  of  the  measurements.  Since  actual  distances  on 
the  target  are  related  to  the  time  delay  by  the  speed  of  light,  the  resolution  is 
approximately  the  speed  of  light  divided  by  the  bandwidth  of  the  measurement. 
It  should  be  noted  that  in  practice  the  frequency  data  is  usually  passed  through  a 
linear  phase  filter  before  the  Inverse  Fourier  Transform  is  taken.  The  filtering  pre¬ 
vents  the  resulting  time  domain  function  from  containing  some  imwanted  ripples 
due  to  Gibb’s  phenomenon. 

For  this  exsunple  a  6”  cube  was  oriented  as  in  Figure  95  and  the  RCS  was 
measured  for  6  =  35.2®,  ^  =  45°  every  10  MHz  while  the  frequency  was  swept 
from  2  GHz  to  18  GHz.  This  data  was  passed  through  a  Kaiser-Bessel  bandpass 
filter  with  parameter  a  =  2.  The  Inverse  Fourier  Transform  was  then  performed. 
The  co-polarized  results  for  a  vertically  polarized  incident  field  are  shown  in  Fig¬ 
ure  104,  and  the  co-polarized  results  for  a  horizontally  polarized  incident  field  are 
given  in  Figure  105.  The  first  pulse  is  due  to  the  scattering  from  the  corner  closest 
to  the  source.  In  both  cases  the  calculations  match  very  closely  with  the  measure¬ 
ments.  The  calculations  match  the  second  pulse  well  for  the  horizontally  polarized 
case,  while  the  measurements  and  calculations  are  significantly  different  for  the 
vertically  polarized  case.  A  first  order  solution  can  only  be  expected  to  match  the 
first  two  pulses  though,  since  the  measured  time  response  includes  higher  order 
terms  which  have  almost  the  same  delays  as  contributions  from  the  other  corners 
making  it  difficult  to  separate  out  the  contributions  from  individual  corners.  The 
contributions  from  each  of  the  illuminated  corners  are  still  difficult  to  identify  due 
to  the  limited  bandwidth  of  the  calculations.  Since  the  bandwidth  of  the  measure¬ 
ments  was  6  GHz  the  resolution  should  be  approximately  0.17  ns.  The  two  way 
time  delays  for  the  illuminated  corners  are  listed  in  Table  1.  Obviously  the  returns 
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Table 


1:  Delay  times  from  different  corners  of  the  6”  cube. 


Corner  No. 

Two  Way  Delay  Time 

(see  Figure  95 

(nsec.) 

,  ,  ' 

1 

0.00 

(ref.) 

2 

0.41 

3 

0.41 

4 

0.83 

5 

0.83 

6 

1.24 

7 

1.24 
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Figure  104:  Impulse  response  for  6”  cube  and  a  6  polarized  source  at  0  =  35.2°, 

<i>  =  45°. 

from  many  of  the  comers  have  essentially  the  same  time  delays  so  there  is  not  a 
separate  pulse  corresponding  to  each  corner. 

5.3.7  Bistatic  Scattering  from  a  2A  plate  illuminated  by  a  fixed  source 

The  scattering  from  a  square  plate  two  wavelengths  on  a  side  is  examined  in 
detail  next.  The  complete  scattering  matrix  (all  four  values  of  <r)  is  found  for  the 
plate  in  the  x-y  plane  with  a  fixed  source  at  —  45°  and  =  0°  (see  Figure  106). 
Circular  pattern  cuts  from  ^  =  0°  to  ^  =  90°  arc  taken  every  30°.  The  results  for 
the  co*potarized  fields  and  the  0°  pattern  cut  are  compared  with  the  previous 
equivalent  currents  and  Method  of  Moment  calculations  in  Figures  107  and  108 
(the  results  for  ^  =  60°  are  given  in  Section  2.7  for  the  previous  comer  diffraction 
coefficients).  The  cross-polarized  fields  are  zero  for  both  a  4  polarized  incident 
field  and  a  ^  pdarized  incident  field.  The  previous  equivalent  currents  and  the 
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Figure  105:  Impulse  response  for  6”  cube  and  a  ^  polarized  source  at  ^  =  35.2°, 

(f>  =  45°. 

new  comer  diffraction  coefficient  agree  so  closely  that  they  essentially  lie  on  top  of 
each  other  in  this  pattern.  This  is  not  unexpected  since,  in  this  pattern  cut,  the 
receiver  is  on  the  Keller  cone  for  the  front  and  back  edges,  which  therefore  give  the 
major  contribution  to  the  scattered  field.  The  new  solution  reduces  to  previous 
solutions  on  the  Keller  cone  (see  Section  4.6),  which  do  not  radiate  any  cross 
polarized  field  for  incidence  normal  to  the  edge.  The  results  for  the  ^  =  30°  pattern 
cut  are  compared  with  the  previous  equivalent  currents  and  Method  of  Moment 
calculations  for  the  co-polarized  fields  in  Figures  109  and  110.  Similarly  the  results 
for  the  cross-polarized  fields  are  given  in  Figures  111  and  112.  The  results  for 
the  ^  =  60°  pattern  cut  are  compared  with  the  previous  equivalent  currents  and 
Method  of  Moment  calculations  for  co-polarized  fields  in  Figures  113  and  114. 
Similarly  the  results  for  the  cross- polarized  fields  are  given  in  Figure  115  and 
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Figure  107:  Co-polarized  RCS  in  the  ^  =  0®  plane  of  a  2\  square  plate  with  a 
polarized  fixed  source  at  —  45®,  =  0®. 


Figure  108:  Co-poiarized  RCS  in  the  ^  =  0®  plane  of  a  2A  square  plate  with  a 
polarized  fixed  source  at  $*  =  45®,  =  0®. 
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Figure  109:  Co-polarized  RCS  in  the  (f>  =  30®  plane  of  a  2A  square  plate  with  a  6* 
polarized  fixed  source  at  =  46®,  <f>*  =0®. 


Figure  HO:  Co-polarized  RCS  in  the  ^  =  30®  plane  of  a  2A  square  plate  with  a 
polarized  fixed  source  at  9*  =  45®,  (f>*  =  0®. 
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Figure  111:  Gross-polarized  RCS  in  the  ^  =  30®  plane  of  a  2 A  square  plate  with 
d*  polarized  fixed  source  at  =  45®,  (f>*  =  0®. 
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Figure  112:  Cross-polarized  RCS  in  the  ^  =  30®  plane  of  a  2A  square  plate  with  a 
polarized  fixed  source  at  =  46®,  <f>^  =  0®. 
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:  113:  Go-polarized  RCS  in  the  ^  =  60°  plane  of  a  2A  square  plate  with  a 
polarized  fixed  source  at  B'  =  45°,  =0°. 


Figure  114:  Co-polarized  RCS  in  the  ^  =  60°  plane  of  a  2 A  square  plate  with  a 
polarized  fixed  source  at  B'  —  45°,  —  0°. 
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Figure  116.  The  results  for  the  ^  =  90®  pattern  cut  are  compared  with  the  previous 
equivalent  currents  and  Method  of  Moment  calculations  for  co>polarised  fields  in 
Figures  117  and  118.  Similarly  the  results  for  the  cross-polariised  fields  are  g^ven  in 
Figures  119  and  120.  Buyukdura  [22]  also  gives  results  using  the  new  solution  for 
the  case  of  a  9  polarised  incident  field  (the  slight  differences  between  his  results  and 
those  given  here  are  suspected  to  be  due  to  a  sign  error  in  his  program).  Overall 
the  new  solution  agrees  well  with  the  Method  of  Moment  calculations  and  does 
not  exhibit  the  false  shadow  boundaries  mentioned  in  Sections  2.6  and  2.7.  The 
previous  equivalent  current  solutions  exhibit  spikes  at  the  false  shadow  boundaries 
which  occur  near  9  —  130®  and  9  =  310®  for  the  ^  =  30®  pattern  cut  and  near 
9  =  240®  and  9  =  300®  for  the  ^  60®  pattern  cut.  In  the  region  from  9  ss  60® 

to  120®  (i.e.  near  the  plane  of  the  plate)  for  both  the  ^  =  30®  and  ^  =  60®  the 
new  solution  and  the  Method  of  Moments  solution  differ  by  more  than  20  dB.  It 
is  suspected  that  most  of  these  differences  are  due  to  the  effects  of  higher  order 
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Figure  116:  Cross-polarized  RCS  in  the  =  60°  plane  of  a  2A  square  plate  will:  a 
<f>'  polarized  fixed  source  at  6*  =  45°,  =  0°. 
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Figure  117:  Co-polarized  RCS  in  the  ^  —  90°  plane  of  a  2A  square  plate  with  a 
polarized  fixed  source  at  —  45°,  =0°. 
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Figure  118:  Co-polarized  RCS  in  the  ^  =  90°  plane  of  a  2 A  square  plate  with  a 
polarized  fixed  source  at  —  45°,  =  0°. 


4^  polarized  fixed  source  at  =  45®,  =  0®. 


terms  (double  and  triple  diffraction,  edge  waves)  which  are  not  included  in  the  new 
solution.  This  may  also  be  the  case  for  in  the  4  =  90°  pattern,  although  this 
also  indicates  that  the  higher  order  terms  for  the  new  solution  and  the  previous 
solution  must  differ  for  points  off  of  the  Keller  cone.  The  results  for  a  pattern  cut 
near  the  plane  of  the  plate  {6  =  39°)  have  been  given  in  Section  4.2. 

6.3.8  Bistatic  scattering  from  a  2A  plate  illuminated  by  a  fixed  source 
in  a  second  position 

In  the  previous  example  the  scattered  field  depended  on  only  one  or  two 
currents  along  each  edge  instead  of  all  three.  This  may  be  easily  seen  from  the 
geometry.  For  a  6  polarized  incident  field  the  incident  electric  field  tangent  to  the 
front  and  back  edges  is  zero  and  the  incident  magnectic  field  tangent  to  the  two 
sides  is  zero.  Therefore  only  Im  ^  edges 

and  If  {Df)  along  the  right  and  left  edges  contribute  to  the  scattered  held.  If  the 

A 

incident  field  is  4  polarized  instead,  the  currents  contributing  from  each  edge  are 
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Figure  121:  2 A  square  plate  in  the  x-y  plane  with  a  fixed  source  located  at 

e*  =  45^  =  45®. 

A 

the  opposite  of  those  for  the  6  polarized  incident  field. 

In  oider  that  all  three  currents  along  each  edge  will  contribute  to  the  scattered 
field  the  source  is  moved  to  a  new  location.  Retaining  the  square  plate  from  the 
above  example  and  moving  the  source  to  9*  -  46®  and  =  45®  (see  Figure  121) 
circular  pattern  cuts  were  once  again  taken  every  30^“  from  ^  =  45®  io  <t>  -  135®. 
The  results  for  the  co-polarised  fields  and  the  4>  ~  ^13°  pattern  cut  are  compared 
with  Method  of  Moment  calculations  in  Figures  122  and  123.  The  croB8-polanz«::cl 

A  A 

fields  are  zero  for  both  a  6  polarized  incident  field  and  a  ^  polarized  incident  field. 
The  results  for  the  (f>  -  75®  pattern  cut  are  compared  with  Method  of  Moment 
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Figure  122;  Co- polarized  RCS  in  the  ^  =  45°  plane  of  a  2 A  square  plate  with  a  (9* 
polarized  fixed  source  at  =  45°,  —  46°. 


Figure  123:  Co-polarized  RCS  in  the  4>  =  45“  plane  of  a  2A  square  plate  with  a 
polarised  fixed  source  at  9'  -  46°,  =  45°. 
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Figure  124:  Co-polarized  RCS  in  the  <i>  =  75°  plane  of  a  2 A  square  plate  with  a 
polarized  fixed  source  at  6^  =  45°,  =  45°. 


calculations  for  the  co-polarized  fields  in  Figures  124  and  125.  Similarly  the  results 
for  the  cross-polarized  fields  are  given  in  Figures  126  and  127.  The  resxilts  for  the 
0  =  105°  pattern  cut  are  compared  with  Method  of  Moment  calculations  for  the 
co-polarized  fields  in  Figures  128  and  129.  Similarly  the  restUts  for  the  cross- 
polarized  fields  are  given  in  Figures  130  and  131.  The  results  for  the  <f>  ~  13t^° 
pattern  cut  are  compared  with  Method  of  Moment  calculations  for  co-polarized 
fields  in  Figures  132  and  133.  Similarly  the  results  for  the  cross-polarized  fields 
are  given  in  Figures  134  and  135.  The  conclusions  which  may  be  drawn  from 
this  example  are  generally  the  same  as  those  of  the  previous  example.  The  new 
solution  and  the  Method  of  Moments  agree  well  overall.  The  region  of  space  where 
the  results  are  not  in  good  ag’eement  is  close  to  the  plane  of  the  plate  as  in  the 
previous  example.  Once  again,  it  is  speculated  that  these  differences  are  due  to 
higher  order  terms  (double  or  triple  diffraction). 


polarized  fixed  source  at  d*  =  45®,  (f>*  =  45®. 


Figure  126:  Cross-polarized  RCS  in  the  ^  =  75®  plane  of  a  2A  square  plate  with  a 
§*  polarized  fixed  source  at  9^  =  45®,  =  45®. 


155 


Figure  127:  Crosa-polarized  RCS  in  the  4>  =  75“  plane  of  a  2A  square  plate  with  a 
polarized  fixed  source  at  9^  —  45®,  =  45®. 


Figure  128:  Co-polarieed  RCS  in  the  4>  =  105®  plane  of  a  2A  square  plate  with  a 
9*  polarized  fixed  source  at  =  45®,  =  45®. 
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Figure  129:  Co-polarized  RCS  in  the  ^  =  105®  pleuie  of  a  2A  square  plate  with  a 
polarized  fixed  source  at  6^  =  45°,  <f>^  =  45®. 


Figure  130:  Cross- polarized  RCS  in  the  <i>  =  105®  plane  of  a  2A  square  plate  with 
a  9'  polarized  fixed  source  at  =  45°,  =  45®. 
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Figure  131:  Cross-polarized  RCS  in  the  =  105®  plane  of  a  2A  square  plate  with 
a  polarized  fixed  source  at  6^  =  45®,  =  45®. 


Figure  132;  Co-polarized  RCS  in  the  <t>  =*  135®  plane  of  a  2A  square  plate  with  a 
9*  polarized  fixed  source  at  9'  =  45®,  -  45®. 
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Figure  134:  Crou-polarited  RCS  in  the  —  135°  plane  of  a  2A  square  plate  with 
a  $*  polarised  fixed  source  at  9*  =  45°,  =  45°. 
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Figure  136:  Cross-polarized  RCS  in  the  4>  =  135°  plane  of  a  2 A  square  plate  with 
a  polarized  fixed  source  at  9*  =  45°,  =  45°. 

5.4  Conclusion 

Several  examples  have  been  given  in  this  chapter  comparing  the  new  corner 
diffraction  solution  with  previous  high  frequency  solutions,  Method  of  Moments 
solutions,  and  measurements.  The  first  set  of  examples  (Sections  5.3.1,  5.3.2,  and 
5.3.3)  illustrate  that  the  new  solution  gives  essentially  the  same  results  as  the 
previous  equivalent  currents  for  patterns  in  the  principle  plane  of  fiat  plates.  The 
next  set  of  examples  for  bistatic  scattering  (Sections  5.3.7  and  5.3.8)  illustrate  that 
the  new  solution  agrees  very  well  with  Method  of  Moments  results  for  a  wide  range 
of  bistatic  an^es.  Farther  most  of  the  differences  between  the  two  are  probably 
due  to  higher  order  terms. 
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CHAPTER  VI 


CONCLUSION 


The  new  equivalent  ctirrents  are  valid  for  perfectly  conducting  structures  with 
straight  edges  and  flat  faces.  The  contribution  to  the  equivalent  edge  currents  due 
to  each  face  of  the  edge  is  given  in  Equations  (3.38)  and  (3.40)  where  the  an¬ 
gles  are  the  edge  fixed  angles  described  in  Figure  1  of  Section  1.5.3.  The  total 
equivalent  edge  currents  for  the  special  case  of  a  flat  plate  are  given  in  Equa¬ 
tions  (3.47)  and  (3.49).  The  equivalent  currents  may  be  used  to  find  the  scattered 
field  in  the  far  zone  when  the  object  is  illuminated  by  a  plane  wave. 

A  new  corner  diffraction  coefficient  for  finite  closed  structures  made  up  of 
perfectly  conducting  flat  plates  has  been  derived  from  the  nev,  equivalent  cur¬ 
rents.  It  has  been  demonstrated  in  Chapter  V  to  agre ;  reasonably  well  with  both 
method  of  moments  calculations  and  measurements,  especially  in  regions  of  space 
where  higher  order  diffracted  fields  are  negligible.  The  contribution  to  the  corner 
diffracted  field  from  one  face  and  one  edge  is  given  in  Equation  (A. 2).  For  the  spe¬ 
cial  case  of  a  fiat  plate  the  contribution  due  to  each  edge  is  given  by  Equation  (A.3). 
The  important  properties  of  the  solution  are  described  in  Chapter  IV.  The  new 
comer  diffraction  coeffident  may  be  used  to  calculate  the  far  zone  scattered  field, 
to  first  order,  from  an  object  illununated  by  a  plane  wave. 

The  new  comer  diffraction  coefficients  and  equivalent  currents  have  some  dis¬ 
advantages  and  some  advantages  over  previous  high  frequency  solutions.  The  main 
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disadvantages  of  the  new  solutions  are  that  they  are  not  valid  in  the  near  zone  and 
some  of  the  terms  must  be  shadowed  in  a  PO  way.  Since  parts  of  the  new  solu¬ 
tions  must  be  shadowed  in  a  PO  way  it  is  unclear  whether  they  can  be  combined 
with  other  UTD  solutions  to  find  the  RCS  of  a  general  three-dimensional  object. 
On  the  other  hand  it  has  been  shown  that  the  new  comer  diffraction  coefficients 
and  the  new  equivalent  currents  give  the  same  results  for  first  order  diffraction. 
The  new  solution  approaches  a  definite  limit  normal  to  the  plate,  reducing  to  the 
well  known  PO  result  in  this  direction.  The  new  solution  also  agrees  closely  with 
Method  of  Moments  solutions  for  bistatic  scattering  from  a  flat  plate. 
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APPENDIX  A 


NEW  CORNER  DIFFRACTION  COEFFICIENTS 


The  new  corner  diflraction  coefficients  are  given  in  a  form  similar  to  previous 
expressions  for  diffraction  coefficients: 
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where  the  plus  or  minus  sign  is  chosen  depending  on  which  endpoint  contribution  is 
being  calculated.  Section  3.5  discusses  choosing  the  correct  sign.  The  expressions 

are  given  by(0-facc  contribution  only) 
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where  the  +  sign  is  associated  with  and  dj’^ 

while  the  -  sign  is  associated  with  the  df*^^,  d^^^,  and  df^  terms.  For  the 
special  case  of  a  Hat  plate  (n  —  2)  the  contribution  from  both  faces  may  be  found 
using 
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APPENDIX  B 


EQUIVALENT  CURRENTS  FOR  EDGE  ON  INCIDENCE 

The  steps  required  to  show  that  the  new  equivalent  currents  axe  finite  for  edge 
on  incidence  (/3'  -♦  0  or  -♦  tt)  are  outlined  in  this  appendix.  Michaeli  [14]  gives 
a  very  brief  outline  of  the  procedure  required  to  find  the  fringe  equivalent  currents 
as  edge  on  incidence  is  approached  in  the  limit. 

Showing  that  the  LPO  and  PO  components  of  the  equivalent  currents  are  finite 
for  edge  on  incidence  is  fairly  trivial,  and  is  only  briefly  described  here.  Using  the 
form  of  the  LPO  components  given  in  Equations  (3.23)  and  (3.24)  and  the  form  of 
the  PO  components  given  in  Equations  (3.29)  and  (3.30),  along  with  the  relations 


Ei 

or  Ejy  sin 

(B.l) 

Hi 

a  iToSin/?' 

(B.2) 

where  Eg  and  Hq  are  the  magnitude  of  the  incident  electric  and  magnetic  fields 
respectively,  it  may  be  easily  shown  that  both  equivalent  currents  approach  finite 
limits  M  -*  0  or  13'  —nr.  It  has  been  assumed  that  the  face  under  consideration 
is  being  illuminated  by  the  incident  field,  otherwise  the  LPO  and  PO  components 
of  the  equivalent  currents  are  identically  zero. 

Finding  the  UTD  components  of  the  equivalent  currents  as  edge  on  incidence 
is  approached  is  more  complicated.  First  an  appoximation  to  the  parameter  a  must 
be  found  for  sm/3'  0.  Using  Equation  (3.45)  and  the  assumption  sin/?^  <  cos/?' 
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results  in 

which  leads  to 


Using  Equations  (B.l),  (B.2),  (B.4),  (B.5),  and  (B.6)  in  Equations  (3.27)  and 
(3.28)  and  then  letting  — >  0  or  letting  if  shows  that  and  are 

finite  for  edge  on  incidence.  The  contribution  from  both  faces  must  be  included  to 
show  that  the  other  component  of  the  electric  equivalent  current  (I^^)  is  finite 
for  edge  on  incidence.  The  contribution  to  the  equivalent  currents  from  the  N-face 
is  found  by  replacing  i  with  —i,  with  ir-0',4'  with  nif  0  with  ir  -  /?,  and 
if>  with  TiTT  —  ^  in  the  expressions  for  the  contribution  to  the  equivalent  currents 
from  the  0-face.  For  the  resulting  expression  for  the  contribution  from  both 
faces  is  given  by 

rVTD  _ 

fcsin/3' 

1  / n(coB  Oo  cot  0’  —  cot  0  cos  ^)  sin  [(?r  ~  Oo)  /n] 
cos  jn)  ~  cos  ((x-  —  Oo)  /n]  sin  ao  *' 
l/n(- cosanCot/9'  -  cot/9cos(nTr  -  <f>))sm  [(tt  -  On)  1^] 
cos  [(uTT  —  ^’ )  /n]  -  cos  [(ir  —  an)  /n]  sin  an 

(B.7) 


a  «  X  -h  j2  In  (sin  0^^ 

(B.3) 

sin  a  « 

1 

(B.4) 

2 

sin^  0' 

1 

2 

sin”  »  0' 

(B.5) 

cos(  ^  )  « 

1 

2 

sm  n  0 

(B.6) 

where  all  of  the  angles  are  measured  with  respect  to  the  edge  fixed  coordinate 
system  associated  with  the  0-face.  For  sin  0  it  is  easily  shown  that  an  «  a© 
and  the  expressions  for  both  an  and  ao  are  the  same  as  the  expression  for  a  given 
in  Equation  (B.3).  Therefore  combining  the  above  with  Equations  (B.2),  (B.4), 


167 


(B.5),  and  (B.6)  in  Equation  (B.7)  one  obtains  an  expression  for  which  is 

valid  for  regions  of  space  where  sin^'  is  small.  After  placing  the  entire  expression 
over  a  common  denominator,  the  limit  taken  as  edge  on  incidence  is  approached 
is  finite. 

All  of  the  equivalent  currents  are  finite  in  the  limit  as  edge  on  incidence  is  ap¬ 
proached  (except  for  tne  forward  scatter  direction).  In  practice  this  means  that  as 
long  as  the  source  is  kept  reasonably  far  away  from  edge  on  incidence  the  computer 
should  be  able  to  take  the  limit  numerically.  The  distance  from  edge  on  incidence 
that  must  be  maintiuned  depends  on  the  precision  of  the  computer  and  which  form 
of  the  equivalent  currents  are  used.  In  this  regard  considerable  errors  may  accur  in 
finding  if  each  face  is  considered  separately,  since  the  contribution  to 
from  each  face  becomes  infinte  as  edge  on  incidence  is.  approached  (although  the 
sum  of  the  contributions  fiom  the  two  faces  remains  finite). 
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APPENDIX  C 

RCS  FOR  NORMAL  INCIDENCE  ON  A  RECTANGULAR  FLAT 

PLATE 


The  RCS  of  a  rectangular  flat  plate  shown  in  Figure  136  is  found  using  the  new 
equivalent  currents.  The  field  radiated  by  the  equivalent  currents  on  an  arbitrarily 
oriented  finite  straight  edge  is  found  first  to  simplify  later  calculations.  The  finite 
straight  edge  is  shown  in  Figure  137  along  with  the  edge  fixed  coordinates  required 
to  define  the  equivalent  currents.  The  currents  are  substituted  into  the  far  zone 
radiation  integral 


(C.1) 


Doing  the  integration  in  closed  form  (only  the  phase  of  El  and  HI  are  functions 
of  z’)  results  in 


sm 


^cos  /3  -  cos  0'^ 


TT  fc  sin /3 sinks' (cos/?  -  cos0') 
|/?  {picBiQ  +  GintZoifto)  - 


(C.2) 


where  the  relations  between  the  equivalent  currents  and  the  functions 
and  Gm  arc  given  in  Equation  (3.53)  and  Equation  (3.54),  while  and  ///q  are 
the  complex  components  of  the  incident  electric  and  magnetic  fields  tangent  to  the 
edge  with  the  phase  reference  taken  to  be  at  the  midpoint  of  the  edge. 
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Z 


Figure  136;  Flat  plate  for  RCS  calculation  with  normal  incidence. 
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Figure  137:  Equivalent  currents  on  a  finite  straight  edge. 

Using  the  above  results  carefully  will  allow  the  RCS  of  the  plate  to  be  found  for 
normal  incidence.  The  edge  fixed  coordinates  for  each  edge  of  the  rectangular  plate 
are  shown  in  Figure  138  along  with  their  relationship  to  the  pattern  coordinates. 
The  relationships  between  the  pattern  coordinates  and  the  edge  fixed  coordinates 
that  will  be  used  later  are 

cos^l  =  sin  dp  sin 
sin  cos  <l>i  ~  —  sin  dp  cos  <i>p 
sin  sin  <f)i  =  cos  dp 

*  COB  ^2  =  —sin  dp  cos  (^p 

8in/32C08^  =  -sindpsin^p 
sin  02  ^  —  cos  dp 

The  expression  for  the  total  scattered  field  may  be  obtained  from  a  superposition 
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Figure  138:  Relationship  between  the  edge  fixed  coordinates  and  the  pattern 

coordinates. 
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of  the  scattered  field  due  to  each  edge.  The  phase  reference  is  placed  at  the  origin 
of  the  pattern  coordinate  system.  Using  the  expression  from  Equation  (C.2)  for 
the  scattered  field  from  each  edge  (in  the  case  of  backscatter,  =  n  —  fi)  results 
in 

1  8in(fct£;co8/3n) 

TT  s  "-1*3  2^  gjjj2  pQg 

{/^n  iGig{^n,<i>n)ElnQ  +  Gin,{/3n>  ^n)2o-fi^ino) 

ina„{ff„,4^)Zo«Lo} 

1  ^  Bin(i:£cos/?n) 

TT  s  ”“^’^2jb8in^^nCOS)9n 
{/?n  (G«(/?n>^n)^'lnO  +  ^n)2o.fftno) 

^nC?m(/3n,<^n)2o.ff^no}  (C.3) 

*2  = 

where  the  exponential  factors  that  have  been  added  shift  the  phase  references  from 
the  center  of  each  edge  to  the  center  of  the  plate  and  component 

of  the  incident  electric  (magnetic)  field  tangent  to  edge  n  at  the  center  of  the  edge. 

Simply  plugging  6p  =  Q  4>n  =  «‘/2,  n  =  1, 2, 3, 4)  in  Equation  (C.3)  will  not 
work  since  some  of  the  terms  become  infinite  here.  However  if  the  limit  is  taken  as 
Bp  approaches  zero,  a  finite  result  may  be  obtained.  The  algebra  required  in  taking 
the  limit  is  fairly  complicated,  so  to  simplify  the  calculations  the  contribution  from 
each  component  of  the  current  (LPO,  UTD,  and  PO)  is  calculated  separately,  and 
then  the  contributions  are  summed  to  give  the  total  scattered  field. 

i»  =  ^  e^TD  _  gPO 


E‘  = 

X 


X 


where 
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(C.6) 

(C.6) 

(C.7) 

(C.8) 


along  with 


Eiio  = 

(C.9) 

miQ  = 

(C.IO) 

- 

(CJl) 

Bizo  = 

(C.12) 

Bho  = 

(C.13) 

^40  = 

(C.14) 

BU^  = 

(C.15) 

=  “4l  h  =  -42 

/^S  =  ^3  =  If- ;  (^1 

^4  =  h  =  '^-<l>2'A<h^^) 

(skipping  the  details)  results  in 

E^PO  ^  ((7.16) 

aIT  s 

As  should  be  expected,  the  same  result  would  be  obtained  using  the  PO  currents 
in  the  surface  integral. 
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Using  the  same  procedure  for  the  UTD  components  of  the  currents,  which  are 


written  in  the  same  form  as  in  Section  4.7, 

ryUTD,^  1  sin ^ sin  I y sin^ /3(1  -  cos <l>)  +  2 cos^ /3 

^te  ^ - .  .  2  a - 1 - 

2  (sin  p  cos  9  —  cos^  p ) 


=  -8in/3co8/3cos|J^ 


2  Y  sin^  /9(1  —  cos  ^)  +  2  cos^  jS 


(C.18) 


=  -;;77 


1  sin^  ^  COB  ■ 


2  (sin^/Jcos^  —  cos^^) 

I  2 


(C.19) 


V  sin^  0{l  —  cos  ^)  +  2  cos^  p 

Plugging  Equations  (C.9)  through  (G.15)  and  Equations  (C.17)  through  {C.19)  in 
Equation  (C.3)  and  taking  the  limit  ta  6p  ~*0  results  in 


^UTD^  +  E*  (C.20) 

IT  S  aIC  $ 

Finally  the  contribution  from  the  PO  components  of  the  currents  is  found. 
Using  the  form  of  the  currents  ^ven  in  Section  4.7 

1  B\n^0Kn<f> 


O  - 

l2  JL  \” _ pl 


2  (sin^  0cos4>~  coe^  /3) 


G£?(fJ.*)  =  0 
G^°{I3A)  =  -G£V.« 


(C.21) 

(C.22) 

(C.23) 


Plugging  Equations  (C.9)  through  (C.15)  and  Equations  (C.21)  through  (C.23) 
into  Elquation  (C.3)  results  in 

=  -  -{wt)- — ^ —  E‘  +  ^(tn  tan^  +  < cot^  ^p)^ — ^ —  £*  (C.24) 

^  s  ^tr  s 

which  is  the  same  as 

Combining  the  results  of  Equations  (C.16),  (C.20),  and  (C.24)  gives  the  total 


scattered  field 


—  ik 


{C.25) 
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or  converting  it  to  RCS 


^co—pol 

^co~pol 


<^x~pol  =  0 


Notice  that  the  fringe  components  (UTD-PO)  cancel  for  backscatter  at  normal 
incidence,  and  the  field  is  determined  by  the  LPO  component  of  the  currents 
(which  predict  the  same  fields  as  the  PO  surface  integral). 
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APPENDIX  D 


AN  EQUIVALENT  EXPRESSION  FOR  THE  EQUIVALENT 

CURRENTS 


The  expressions  for  the  equivalent  currents  in  Equations  (4.16)  and  (4.18)  are 
derived  here  from  previous  expressions  in  Buyukdura  [22]  and  Michaeli  [13].  Since 
the  expressions  are  quite  complicated  each  of  the  three  components  (LPO,  UTD, 
and  PO)  will  be  derived  separately. 

Starting  with  the  LPO  components  of  the  currents  from  Buyukdura  [22], 


jLPO  _ 
j^LPO 


210 

jk  ^*8in/3' 
2Zq  i  sin  ^ 
jk 


2 

C  +  — (cot  /3  cos  ^  -f  cot  13*  cos 
jk 

C 


(D.1) 

(D.2) 


where 


C  = 


sin  0  cos  (j>  +  sin  cos  (f>' 


{D.Z) 


[(cos/?  -  co8/3')2  -f  (sin/? cos  ^  +  sin^'  cos^')^j 
Farther  separating  the  equivalent  electric  current  into  and  components 
results  in 


tLPO  ^  2Kop<rin^' 


jk  “‘sin/S' 
combining  with 


/?  cos  ^  +  sin  /?♦  cos  ^) 


(D.4) 


coa^l  =  sin^isin^cos^  +  cos^l  cos/? 


(D.5) 
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co6^1  =  sin^i  sin^^cos^^  4- CO8  01  cos;9^ 

^  „  cos  /9  —  cos  0' 

cot  6l  =  -; - - - - - 

sin  0  cos  ^  4-  sin  cos  (j>' 

results  in 

jLPO  _  2Yq ^ism<f>'  eiR 9 1 _ 1 _ 

^  jk  ^  sin/S'  cGsV*!  +  co8^»j 

Using  a  trigonr“-  ;tric  identity  from  Section  3.4  results  in 

tLPO  _  i^O  sin  sin  ^1 
*  2k  ^sin/J'sinV'j 

I  |eot  _  „,t  (Z±JiLz£l\ 


(D.6) 

(D.7) 


(D.8) 


+  (V>1  +  V’i) 

4 


The  expression  here  is  the  same  as  that  given  in  Section  4.7  except  for  the  step 
function  which  is  implied  in  this  expression  since  Buyukdura  assumes  0-face  in¬ 
cidence.  Using  a  similar  procedure  the  expressions  in  Section  4.7  for  and 
ill  may  be  derived. 

Next  the  expression  for  is  developed.  Starling  with  the  expression  for 
derived  by  Buyukdura  (22). 
rVTD  2j‘  i  sin/3'co6^ 


”*  k  *sin^'(Bin/3co8^  4- sin/S^  cos^) -t- (cos^  -  cos^')co8/9^ 

X  jeot  0*  (sin  ^  cos  4  (cos /(?  —  cob^)  cot  5*)  -  sin  cot  (3  cos^j 

sin/J'cos^')  -  (cos/3  -  coej9')coaj3' 

(although  the  terms  are  not  separated  in  this  way  by  Buyukdura  [22],  it  cair  be 
shown  that  the  part  of  1  given  here  corresponds  to  ) 


Using 


cos  J2  -  cos  0' 

cos  V’S  =  sin  ^  sin  0  cos  4>  +  cos  62  cos  0 


cosV’2  =  sin 62 cos 4  €06^2 cos/?* 
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(D.IO) 

(D.ll) 

(D.12) 


results  in 


rUTD 


‘m 


2  I  sin/fl'cos^  .  I 
jk  *  cos  tj}2  +  cos  V  1  — 


where 


cot/3'  r  2af] 

*"»  =  -  cos  ^  ]  - 


if}'2  Vl-C08V’2 

sin  /3'  cos  ^  cos  ^ 


sin /S'  1  J  sin/S 

Using  the  trigonometric  identity  as  before  results  in 

J.I 

‘m  ~  riJ  -  ./.»  . 


(D.14) 


jVTD  _  jg,  irinff'cos^-  . 

sin  ^  sin  ^2 

X I  |coi  j  -  col  ^ 


+ 


cot 


-(V>2  +  V>2)j  _  ^^r  +  (V>2-<-V>2)y 


(D.15) 


which  is  the  some  as  the  expression  pven  in  Section  4.7.  The  results  for  and 
^jVTD  obtained  in  a  similar  manner. 

Finally  the  expression  for  of  Section  4.7  is  found.  Starting  with  the 
expression  given  by  Michaeli  [14| 

PO  2^0  rrtrr/.^  il\ 


M 


jk  ^  ^  sin/?  sin /3' (cos  cos  a) 


sin  /3  cos  4>  (cos  8  —  cos  /3')  cos  /?' 
cos  a  =  — : — T, —  + -  - 


sin  0'  '  sin^/?' 

Substituting  the  expression  for  cos  o  into  the  expression  for  results  in 


(D.16) 

(D.17) 


M 


=  ^g;c/(r .. 

jk  ^  ^  »n0 


sio'^/3'coe^ -f  wn/S'Mn/3co8  0  4- (cos/?  -  co6^')co5;3'  ^  ^ 

Using  ^2*  and  02  defined  previously. 

1 


jk  amp  c«KiV2  +  codv^ 


(X).19) 
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Using  the  trigonometric  identity  from  Section  3.4  as  before  results  in 


^-i^2+i>2)^  _  +  (V>2+V>2)y  I 


izlj 


+ 


cot 


(D.20) 


which  is  the  same  as  the  expression  given  in  Section  4.7.  The  apparent  difference 
due  to  the  replacement  of  the  5*  function  by  the  step  function  f^(ir  -  is  only  a 
technicality  since  the  expression  given  here  gives  the  contribution  from  boih  faces, 
but  is  valid  only  for  O-face  incidence  while  the  expi^ssion  in  Section  4.7  gives  the 
contribution  from  both  faces  for  either  O-face  or  N-face  incidence.  The  other  two 
terms  and  may  be  found  similarly. 
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